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ABSTRACT 

■^r I Context. Because of their relatively simple morphology, "bubble" H II regions have been instrumental to our understand- 

r^ , ing of star formation triggered by Hll regions. With the far-infrared (FIR) spectral coverage of the Herschel satellite, 

^^ ■ we can access the wavelengths where these regions emit the majority of their energy through their dust emission. 

r^ ' Aims. We wish to learn about the dust temperature distribution in and surrounding bubble H II regions and to calculate 

^ ll the mass and column density of regions of interest, in order to better understand ongoing star formation. Additionally, 

I , we wish to determine whether and how the spectral index of the dust opacity, /3, varies with dust temperature. Any 

O ■ such relationship would imply that dust properties vary with environment. 

.fj ' Methods. Using aperture photometry and fits to the spectral energy distribution (SED), we determine the average 

^2 ' temperature, /3-value, and mass for regions of interest within eight bubble Hll regions. Additionally, we compute maps 

^ of the dust temperature and column density. 

Results. At Herschel wavelengths (70 /^m to 500 fim), the emission associated with HlI regions is dominated by the 
cool dust in their photodissociation regions (PDRs). We find average dust temperatures of 26 K along the PDRs, with 
little variation between the HlI regions in the sample, while local filaments and infrared dark clouds average 19 K and 
15 K respectively. Higher temperatures lead to higher values of the Jeans mass, which may affect future star formation. 
»vj ' The mass of the material in the PDR, collected through the expansion of the Hll region, is between ~ 300 Mq and 

~ 10, 000 Mq for the HlI regions studied here. These masses are in rough agreement with the expected masses swept 
up during the expansion of the H II regions. Approximately 20% of the total FIR emission is from the direction of the 
bubble central regions. This suggests that we are detecting emission from the "near-side" and "far-side" PDRs along the 
CO ■ line of sight and that bubbles are three-dimensional structures. We find only weak support for a relationship between 

^3 ' dust temperature and /?, of a form similar to that caused by noise and calibration uncertainties alone. 
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1. Introduction ([Elmegreen fc Lada]ll977D . It has been shown in simulations 
C^ that this layer may contain several thousand solar masses 
- - Hii regions that appear as a rmg at mfrared (IR) waye- (|Hosokawa fc Inutsukal l2006l) . a fa ct that has been con- 
lengths, or a bubble seen m projection, have been the firmed for ind ividual Hii regions (iDeh arvcng et all [20031: 
focus of numerous studies of triggered star formation be- IZavagno et aLl[200l [20071: IPoSareset a l. 2009). 
cause of their relatively simple morphology. Hii regions 

expand as they age due to the pressure difference be- The bubble morphology is common for Galactic Hii 
twee n the ionized gas and the surrounding neutral medium regions. IChurchwell et al.| (|2006|, 12007D compiled a cata- 
(see iDvson fc Wilfe[^ [Tool . During this expansion, a log of -^ 600 bubbles detected at 8.0 /xm in the Spitzer 
layer of collected neutral material can form on the bor- Galactic Leg acy Infrared Mid-Plane Survey Extraordinaire 
der of the Hii region; within this collected layer new stars (GLIMPSE; [Beniamin et^ I2003D. Nearlv all identi- 
may form. This process is known as "collect and collapse" fied bubbles enclose Hii regions (|Deharveng et al.| |2010|; 
iBania et al.ll2010l : lAnderson et al.ll2011l ). and furthermore. 



* Herschel is an ESA space observatory with science instru- nearly half of all H II regions have a bubble morphology 
ments provided by European- led Principal Investigator consor- ([Anderson et al.l [20 111) . Most of these bubbles have sur- 
tia and with important participation from NASA. rounding material that appears to have been collected 
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through their expansion (jDeharveng et al]l2010t ). Knowing 
the temperature of this material is necessary to better esti- 
mate the mass and column density in the collected layers, 
which in turn is necessary to determine the efficiency of 
triggered star formation. 

Although it contains just ~ 1% of the mass, dust plays 
a significant role in the energetics of Hii regions. Dust acts 
as a coolant for Hii regions - it absorbs high energy pho - 
tons and re-emit s in th e IR. IWood fc Churchwell (|l989b[ ) 
and iKurtz et al.l (|1994f ) estimate that for ultra compact 
(UC) Hii regions, dust absorbs between 42% and 99% of 
the ionizing photons. The absorption by dust leads to a 
slower expansion rate and a smaller physical size, stalling 
the expansion of an Hii region earlier than it would other - 
wise. This phenomenon has been stud ied bv lMathisI (Il97ll). 
iPetrosian eFaLl (|1972t) , ISpitzeiJ (|1978[) . and by 
(1200411 . 



The derived dust properties tell us about the grain pop- 
ulation, which in turn tells us about the physical conditions 
of the environment. Aside from the temperature, we may 
also learn the wavelength dependence of the dust opacity, 
K^. This re lationship can be modeled as a power law (cf. 
lHildebrandl [r983): 



Ky = K(3{v/V()Y ■ 



(1) 



Arthur et al.l 



Observations of H ii regions from mid-IR to mm- 
wavelengths trace the re-radiated energy from dust. Such 
observations can be used to derive the column density and 
mass distributions. Using observations of dust to estimate 
the total mass of gas and dust has the advantage of being 
applicable over a large range of column densities. Optically 
thin far-IR (FIR) to mm-wavelength observations of dust 
can be used to trace the mass distribution in high-density 
environments where low-density gas tracers such as CO 
would freeze out onto dust grains, and also in low-density 
environments where high-density gas tracers would not be 
detectable. 

Until recently, the dust temperatures of H ii regions were 
poorly constrained. Previous authors have used data at 12, 
25, 60, and 100 ^m from the Infrared Astronomical Satellite 
(IRAS) to derive flux ratios (colors), and thus infer a dust 
temperature. These studies found that for nearly all H ii re- 
gions the flux at 100 /im is greater than the flux at 60 /^m, 
implying dust temperatures ^ 30K. This phenomenon oc- 
curs for giant H ii regions (those with Lyman continuum 
photo n emission rates iViyc > 10^° s~^) (Conti & Crowthefl 
I2004D . optically visible Hii regions (fchan fc Fich 1995^ 
and ultra compact (UC) H ii re gions fiWood fc Churchwell 
Il989bt ICrowther fc Contill2003t ) One of the warmer H ii re- 
gions known is Ml 7, for which [Povich et all (|2007fl found 
that the 60 /im flux is greater than the 100 /xm flux. With 
one data point on the Rayleigh- Jeans part of the spec- 
tral energy distribution (SED), they were able to better 
constrain the dust temperature. They found a dust tem- 
perature in the photodissociation region (PDR) of M17 of 
~ 40K and a dust temperature of the "interior" in the 
direction of the ionized gas of ^ 100 K. 

With the advent of the Herschel Space Observatory 
(jPilbratt et al. 2010), we have access to the crucial FIR 
regime long-ward of the IRAS bands, at high angular res- 
olution. At the longest IRAS wavelength band of 100 /im, 
the in-scan resolution was ~ 3' and the out-of-scan reso- 
lution was ~ 4'. The photometric bands of Herschel span 
70 /im to 500 /xm, at resolutions from 6" to 37". Thus, while 
only very rarely were H ii regions observed to have an IRAS 
data point on the Rayleigh- Jeans side of the SED, Herschel 
can provide a well-sampled SED with multiple data points 
in the Rayleigh- Jeans side. Furthermore, the resolution of 
Herschel allows us to determine dust temperature varia- 
tions within H ii regions, which was largely not possible 
with IRAS. 



where kq is the opacity at frequency vq and j3 is the spec- 
tral index of the dust opacity. This relationship does not 
hold a shorter wavelengths in the mid-infrared (MIR) where 
the dust emission profi le is more complicated (see model in 
ICompiegneet al.ll20lH) . 

Many authors have empirically found support for an 
inverse relationship between the dust temperature and (3. 
For example, such a relationship was found in PRONAOS 
observations of Ga lactic cirrus and star-forming regions 
(jDupac et al.ll2003[ ). ARCHEOPS sub-mm point sources 
(including numerous star formation regions) (jDesert et al.l 
200^, BOOMERANG observations of Galactic 



cirrus 



(Ve neziani et al.l l2010t) . Herschel Hi-Gal (Molinari et al. 



2010) observations of the Galactic plane (Paradi s et al, 
2010 ), PLANCK observatio n s of ne arby molecular clouds 
f|Planck Collaboration et al .' l2011ajl, PLANCK observa- 



tions of cold cores (iPlan ck Coll aboration et al. 2011bD. 
BLA ST observations near the Vela complex (jMartin et al.l 
1201 ll ). and ind ividual Hii region environments observe d 
with Herschel (jAnderson et al.l 120101: iRodon et aLll2010[ ). 
It has also been fou nd in laboratory work of dust similar 
to th at of the ISM (jMenneUa et al.l(l998l: ICoupeaud et al.l 
[201l|). Such a relationship, however, can be falsely pro- 
duced by measur ement noise and lin e o f sight temper- 
ature variations (IShettv et al.l l2009al lbl: iJuvela fc YsardI 
I2012D . iMalinen et al.l ()2010[ ) find that a false relationship 
is caused in the vicinity of a heating source. Classically, 
we would expect that /3 i s between 1.0 and 2.0 (see 
iTielens fc Allamandolalll987l and references therein) - for 
a blackbody, /? is equal to zero. The (i — Td relationship 
has been proposed to arise naturall y from the disord ered 
structure of amorphous dust grains ()Meny et al.ll2007f ). 

Us ing Herschel sc i ence demonstration phase (SDP) 
data, lAnderson et al.l ()2010t) . derived dust temperatures 
and dust /3- values for the Galactic Hii region RCW120. 
They found temperatures of the dust associated with the 
Hii region from 20 to 30 K, with colder nearby patches 
of 10 K associated with in frared dark c l ouds (IRDCs). A 
similar result was found bv lRodon et al.l (|2010[ ) for Shl04, 
again using Herschel SDP data. The results of both papers 
were consistent with a relationship between /3 and Td, with 
the coldest regions having /3- values near 3.0. In the present 
work, we extend these analyses to a sample of eight "bub- 
ble" H II regions observed by Herschel. 



2. HII Region Sample 

Our sample includes eight bubble H ii regions: Sh 104, W5- 
E, Sh241, RCW71, RCW79, RCW82, G332. 5-0.1, and 
RCW120. With the exception of W5-E and G332.5-0.1, 
all H II regions in ou r sample were list e d as c ollect and col- 
lapse candidates in iDeharveng et al.l (|2005t ). All but one 
region, G332.5— 0.1, have been identified from their opti- 
cal emission and thus extinction and the amount of in- 
tervening dust along the line of sight is low. The basic 
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parameters for the regions are given in Table 1, which 
hsts the right ascension, dechnation, Galactic longitude, 
and Galactic latitude of the approximate center position, 
the approximate angular diameter, the assumed distance, 
and the physical diameter. Our sample targets span a 
range of angular diameters from 3' to 34', distances from 
1.3 kpc to 4.7 kpc, and physical diameters from 2pc to 
20 pc. Three-color Herschel images composed of data from 
500 /zm, 250 //m, and 100 /zm respectively in the red, green, 
and blue channels are shown in Figure [T] these data are 
discussed in 521 The contours in Figure [1] show the ra- 
dio cont inuum emission from the NRAO VLA Sky Survey 
(NVSS; ICondon et aJlllQQSD. Syd ney University Molonglo 
Sky Survev CSUMS S:lBock et al1ll9 99'). or the Green Bank 
6 cm survey (GB6: iGregorv et al]|l996i ). The radio contin- 
uum emission traces the ionized gas of the Hii regions. 
These contours are only meant to show the strongest radio 
continuum components. Due to the sensitivity of these sur- 
veys, much of the diffuse emission is not detected, including 
in some cases emission from the interior of the bubbles. 

Throughout we refer to the direction of the central re- 
gion of the bubbles as the "interior," and the collected ma- 
terial delineating the bubble structure as the "PDR." For 
all regions in our sample, the interior is spatially coinci- 
dent with the ionized gas. We use the term "IRDC" for 
any filament that is detected in absorption at 8.0 /xm. 

2.1. Shl04 



Shl04 (|SharDlessl[l959D is a -- 7'-diameter Hii region. It 
has a shell of material located along th e PDR that was 
detect ed in CO (2 — > 1) observations by iDeharveng et al.l 
(|2003f) . They estimate that the total mass in the shell is 
~ 6, 000 Mq. CS observations of the same field show four 
regularly-spaced dense molecular clumps on the border of 
Shl04. One clump, of mass ~ 670 M©, is spatially coinci- 
dent with an Hii region detected in the NVSS. This region 
is seen on the eastern PDR of Sh 104 in Figure [TJ 

iRusseill (|2003f ) lists a kinematic distance for Shl04 of 
4.5lo'8kpc and a spec troscopic, or spe c tro-ph otometric, 
distance, of 3.1ib0.9kpc. lDeharveng et al.l(|2003D used4kpc 
in their analysis, which we adopt for the present work . 
Sh 104 is exci ted by an 06V star (jCrampton et al.lll978t 
lLahullalll985[ ). Given the angular diameter, the physical 
diameter is ~ Spc. 

Sh 104 has strong 100 /xm emission in its interior region 
('see l5.4|) . Also, it is the most "complete" of the Hii regions 
in the sample in that the PDR is continuous and surrounds 
the entire region. These two phenomina may be related as 
the completeness of the PDR may indicate that radiation 
pressure and stellar winds are not particularly strong. In 
this scenario, the increased emission at 100 /xm is a sign 
that the radiation is trapped in the bubble interior. There 
are numerous cold filaments that begin on the border of 
Sh 104 and lead radially away. 

2.2. W5-E 



W5-E ()Westerhoutlll958D is a '- 34'-diameter outer Galaxy 
H II region that has been the focus of numerous studies. It 
is part of the W3-W4-W5 complex; W4 and W3 are to the 
west of W5. Our observations are of W5-East, one part of 
the W5 region; the other part is W5-West, which is a second 



bubble that appears to be interacting with W5-East. Here 
we refer to "W5-East" as "W5-E." 

Here we adopt a distance of 2.0 kpc for W5-E. This 
distance is the same as the mase r parallax distance as - 
signed to W3(0H), 2.04±0.07kpc (|Hachisuka et al.l[2006h . 
and agrees well wit h estimates from spectroscopic paral- 
lax measurements ([Becker fc FenkartI 119711: iMoffatI Il972l: 
Massev et al.l Il995fl . The main exciting star of ^y5-E is 
HD 18326, which is of spectral type 07V (|Walbornlll97l . 

The Herschel image in Figure [T] shows a nearly complete 
bubble, open to the south. This southern opening is likely 
due to a density gradient in the ambien t medium, as can 
be inf erred from the CO observations in lHever fc TerebevI 
(Il998l ). There are also numerous "bright- rimmed clouds" 
(BRCs) and "elephant trunks" pointing toward the bub- 
ble interior. These features were seen in Spitzer observa- 
tions of the region (jKoenig et al.ll2008l ). and are prime sites 
for investigating triggered star form ation. Our obse r vation s 
contain three BRCs cataloged by ISugitani et al.l (|1991[ ): 
BRC12, BRC13, and BRC14. 

Point sou r ces d etected by Spitzer were analyzed by 
iKoenig et al.l ()2008D . who found two distinct populations 
of protostars; a population of Class II young stellar objects 
(YSOs) toward the interior of W5-E and a population of 
Class I YOSs in the surrounding molecular material. They 
hypothesize that these populations are separated by age 
and that triggering may have caused the creation of the 
Class I population. 

A detailed study of star formation in W5-E using 
these same Hers c hel da ta is given in a companion paper, 
IDeharveng et al.l (|2012f ). 

2.3. Sh241 



Sh241 ()Sharplesslll959t ) is a ~ 3'-diameter Hii region. The 
region itself shows only faint Ha emission (Pomares et al., 
2012, in prep.), although there is a dense adjacent molecu- 
lar core that has been the focus of numerous studies in CS 
and HON (Plume et al.l ll992t iPirogovl [l999l: IShirlev et all 
2003; Wu et al. 2010). It has also been mapped at 350 /itm 
(M ueller et al.ll200^ . a molecular outfiow has been detected 
( . Wu et al.l I19 99I). and an H2O maser has b een detected 



Cesaroni et al. 1988- iHenning et all I1992D . iMoffat et ail 



19791 ) find that the spectroscopic distance to Sh241 is 



4.7 kpc. 

Sh241 is the least complete of the bubble Hii regions in 
our sample. Figure [T]shows that Sh 241 is open to the south 
and defined in the north by material that emits strongly in 
the FIR. There are numerous condensations in the field, 
especially to the west and south, and there is a separate 
Hii region to the west detected in the NVSS. We note that 
there is faint emission seen to the south in Figure [T]that 
may be part of a secondary PDR of Sh241. It is unclear if 
this is truly a secondary PDR however and the reported 3' 
size refers only to the more compact emission seen in at the 
center of the Herschel data. 



2.4. RCW71 



ROW 71 (jRodgerset al.lll960[ ) is a ~ 3'-diameter Hii re- 
gion. Compared to the other H 11 regions in the sample, little 
is known about RCW71. The Herschel images in Figure |T] 
show numerous small condensations to the south and west 
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Table 1: Basic properties of the observed Hii regions 



Name 


RA (J2000) 


Dec. (J2000) 


; 


b 


d 


Ang 


. Diam. 


Phys. Diam. 




hh:inni:ss 


hh:inin:ss 


deg. 


deg. 


kpc 


arcinm. 


pc 


Shl04 


20:17:42 


+36:45:26 


74.761 


+0.618 


4.0 




7 


8 


W5-E 


02:59:41 


+60:32:01 


138.073 


+ 1.488 


2.0 




34 


20 


Sh241 


06:03:58 


+30:15:25 


180.872 


+4.110 


4.7 




3 


6 


RCW71 


12:50:21 


-61:34:57 


302.803 


+ 1.289 


2.1 




3 


2 


RCW79 


13:40:17 


-61:44:00 


308.681 


+0.592 


4.0 




10 


12 


ROW 82 


13:59:29 


-61:23:40 


310.984 


+0.409 


3.4 




6 


6 


G332. 5-0.1 


16:16:59 


-50:48:14 


332.523 


-0.136 


4.2 




4 


5 


ROW 120 


17:12:24 


-38:27:44 


348.252 


+0.474 


1.3 




8 


3 




Fig. 1: Three-color Herschel images for the eight observed regions. For each source the 500 /im data are shown in red, the 
250 /im data in green, and the 100 /ini data in blue. The images are oriented in RA and Dec. such that north is up and 
east is left. The white scale bar in the upper left corner is 5' in length. The coldest regions in each field appear red while 
the warmest appear blue. Contours show radio continuum emission from the NVSS (for Shl04, Sh241, and RCW120), 
SUMSS (for RCW79, RCW82, and G332) or GB6 (for W5-E and RCW71). 



of RCW 71 that are bright at SPIRE wavelengths. The CO velocity of these condensations is similar to that of RCW 71 
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Fig. 1: - continued. 



and are therefore likely associated (M. Pomares, private 
communication). The PDR of RCW71 seen in Figure [T] is 
filamentary. The spectroscopic distance of the exciting star 
of RCW71, HD311999, is 2.11 kp c while its ne ar /far kine- 
matic distances are 3.05/6.15 kpc ('Rus seill[2003f ) . The spec- 
troscopic distance is in better agreement with dist ances to 
nearby young stellar clusters ((Russeil et al.lfl998r) and we 
adopt the spectroscopic distance of 2.1 kpc for the present 
work. This gives a physical diameter for RCW 71 of '^ 2 pc. 



mass of the most m assive condensation is ^ 1, 000 Mq. 
IZavagno et al.l (|2006[ ) detected several Class I YSOs in the 
most massive condensations. The bright compact source 
seen to the south-east in Figure [T] is a separate H ii re- 
gion detected in radio continuum emission with SUMSS. 
Figure [1] shows that RCW 79 is open to the northwest. To 
the south, there is what appears to be a second ionization 
front. 



2.5. RCW 79 



RCW 79 (jRodeers et al.lll960[ ) is a ~ 12'-diam eter Hii re- 
gion . Using 1.2 mm continuum observations. Z avagno et al.l 
(|200 6') found a fragmented layer of neutral material along 
the PDR of RCW 79 with a total mass of ~ 2, 000 Mq . The 



iRusseill ()2003[ ) find a kinematic distance of 4.8]tiokpc 
for RCW 79 and a spectroscopic distance of 4. 0±0. 6 kpc. We 
adopt 4 kpc for t he present work, whic h leads to a physical 
diameter of 12 pc. iMartins et"all ()2010D found that RCW 79 
is ionized by a cluster of a dozen O stars. 
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2.6. RCW82 



RC W82 (iRodgers et al 



■ :. -- 19601 ) is a 6'-diaineter Hii re- 
gion. iPomares et al.l (J2009) studied the molecular emission 
traced with CO isolopologues surrounding RCW82 and 
found a fragmented shell of emission with a total mass 
of ~ 10, 000 Mq. The most massive condensations have 
masses of ~ 2 500 M p) . The YSO population identified by 
iPomares et al.l (J2009t) is not evenly distributed, but is con- 
centrated on the border of RCW82, indicating that trig- 
gered star formation may have lead to their formation. 

RCW82 appears along an IRDC filame nt detected on 
both s ides of the H ii region. The CO data in IPomares et al.l 
(|2009tl show that this filament is at the same velocity as 
RCW82, and that the filament and Hii region are there- 
fore associated. There are numerous other filaments seen in 
Figure [1] that lead radially away from RCW82; most are 
not seen in absorption at 8.0 /im. 

The r e is so me uncertainty as to the distance to RCW 82. 
iRusseill (|2003( ) lists a kinematic distance of 4 3±0.7kpc and 
a spec troscoptic distance of 2.9 ± 0.9 kpc. IPomares et al.l 
(|2009tl use a distance of 3.4 ± 0.9 kpc, which we adopt for 
the present work. Given the angular diameter and this as- 
sumed distance, the physical diameter of RCW 82 is 6pc. 
iMartins et"an (|2010f ) find that RCW 82 is ionized by two 
09-B2V/III stars. 



2.7. G332.5-0.1 

G332.5— 0.1 is a 3'-diameter Hii region. It is located along 
a prominent IRDC running east-west, detected in emission 
at longer wavelengths with Herschel, that is associated with 
G332.5— 0.1 in velocity (Deharveng et al., 2012b, in prep.). 
This filament, seen in Figure [I] has numerous sources 
within it detected by Herschel. G332.5— 0.1 is unique in our 
sample in that hot dust at 100 fim is faint in its interior. 

On the northern border of G332.5— 0.1 there is an as- 
sociated UC Hii region that has been the focus of nu- 
merous studies_as_£ajrtof the Red MSX Source survey 
fRMS:' Urauhart et al.l[2008[) . Just off the western border of 
G332.5— 0.1 there is another region of extended radio con- 
tinuum emission detected at 843 MHz with SUMSS and a lso 
at 24 fim with the MIPSGAL survey (|Carev et al.ll2009D - 
it is likely a distinct compact Hii region associated with 
G332.5-0.1. 

The recombinati on line velocity o f G3 32.5— 0.1, 
—50 kms^^ from I Caswell fc HavnesI (|l987l ) places 
G332.5— 0.1 at a kinematic distance of 3.7 kpc. This 
assumes the near kinematic distance, which seems likely 
given its association with IRDCs (IRDCs at the far 
distance would be difficult t o detect due to a lack of 
background behind the cloud). JRusseil et all ()2005D have a 
slightly revised distance of 4. 2 kpc for G332.5— 0.1 based on 
nearby Hii regions with a similar velocity; we use 4.2 kpc 
here which gives a physical diameter of ~ 5 pc. 

2.8. RCW120 



RCW 120 (jRodgers et al.l ll960D is an 8'-diameter H ii re- 
gion that has been the focus of numerous recent stud- 
ies. It has a massive, fragmented lay er of neutral materia l 
seen along its PDR traced at mm (jZavaeno et al.l 120071 ) 
and sub- mm (jDeharveng et al.l I2009D wavelengths. Of th e 
eight mm-condensations located by IZavagno et al.l (J2007l) , 



five are found on the PDR, indicating that this material 
has bee n collected duri ng the expansion of the H ii region. 
Deharv eng et al.l (|2009[ ) estimated that the mass of the col- 
lected layer is -- 2, 000 Mq. 

RCW 120 i s among the nearest H ii regions to the Sun. 
IRusseill ()2003() lists a kinematic distance of 1.8^0^7 kpc 
and a stellar distance of 1.3 ± 0.4 kpc for RCW 120; we 
adopt 1.3 kpc for the present work. Given the angular size, 
the physical size of RCW 120 is ^ 3pc. A single star is 
i onizin g RCW 120: CD-38 "11636. iGeorgelin fc GeorgelinI 
(|1970[ ) found that its spect ral type is 08 f rom spectro- 
scopic measurements while (jCramptonl 1197 1") find that it 
is an 09. More recently , using near-IR integral field spec- 
troscopy, IMartins et"aI1 (|2010l ) found that RCW 120 is ion- 
ized by a 06-8V/III star. 

There is significant active star formation in the 



Zavagno et all (I2007D . 



surroundings o f RC W 120. ^, _ , 

iDeharveng et al.l (|2009D . and IZavagno et al.l (|2010D all 
located YSOs in the field. The most massive condensation 
is known as "Condensation 1" and is loca ted to the south- 
west i n Figure [H Within Condensation 1. [Deharveng et al.l 
(|2009l ) found an evenly spaced chain of eleven Class I or flat 
spectrum sources parallel to the IF. The authors suggest 
that it is an example of the Jeans gravitati onal instability. 
Using the same Herschel data shown here, 'Zavagn o et al.l 
(201(|) detected a massive YSO towards Condensation 1 
with a stellar mass of 8-10 M©. This c ondensation was 
detec t ed previously a t TOjxm by Spitzer (jDeharveng et al.l 
I2009D . IZavagno et al] (|2010D suggest that this is the first 
detection of a massive Class object formed by the collect 
and collapse process on the border of an H 11 region. 

There are also numerous IRDCs seen in the periph- 
ery of RCW 120, many with point sources observed in 
the mid-infrared suggestive of active star formation. The 
largest IRDC, G348.40-H00. 47, is observed to th e north- 
east of RCW 120 in FigurefTl lJackson et al.l ()2008f ) measure 
a CS velocity of -6.8 kms"! for G348.40-h00.47, which is 
close to the recombination line velocity of the ioniz ed gas 
in RCW120, -12 kms'^ (|Caswell fc HavnesI [19871) . The 
IRDC and the H 11 region are therefore likely associated. We 
assume the smaller IRDCs to the north and west are also 
at the same distance as RCW 120. There are numerous cold 
filaments to the south of RCW 120 observed by Herschel, 
some of which are also seen in absorption at 8.0 /xm. 

3. Observations and Data 

In addition to the Herschel data that is the main focus 
of this work, we utilize auxilliary data from Herschel Hi- 
Gal, Spitzer GLIMPSE, Spitzer MIPSGAL, and APEX- 
LABOCA ATLASGAL. These auxfllary data are available 
for four of our objects: RCW 79, RCW 82, G332.5-0.1, and 
RCW 120. 

There are numerous point sources detected at all wave- 
lengths for the H 11 regions in our sample. We make no at- 
tempt to remove them from the data. This is primarily be- 
cause we cannot with confidence separate the sources from 
the local background at all wavelengths, while ensuring that 
we are locating the same source at all wavelengths. For ex- 
ample, an embedded source in a small IRDC can easily be 
identified at MIR wavelengths, but at the longest Herschel 
wavelengths with lower angular resolution, its emission ex- 
tends over the entirety of the IRDC. If we were to remove 
this source, we would likely be removing a large fraction of 
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the flux from the IRDC as weU, preferentially for data with 
lower angular resolution. 

3.1. Herschel HOBYS and "Evolution of Interstellar Dust" 

Our sources were observed by the Herschel Space 
Observa tory with the PAC S (|Poglitsch et all |2010D and 
SPIRE (iGrifRnet all 120101) instrumentf] as part of the 
HOBY S (iMotte et al.lfeoiOf) and "Evolution of Interstellar 
Dust" (jAbergel et al.ll2010l ) guaranteed time key programs. 
Data were taken in five wavelength bands from 100 fim to 
500 urn: 100 fim and 160 /xm for PACS at FWHM resolu- 
tions of 6.7" and 11", and 250 fim, 350 /zm, and 500 /xm for 
SPIRE at FWHM resolutions of 18", 25", and 37". For both 
instruments, the observations were composed of two orthog- 
onal scans. For all observations with the PACS instrument, 
the scan-speed was 20" per second while for SPIRE it was 
30" per second. The observational parameters are summa- 
rized in Table 2, which lists for each source the field size, the 
total integration time, the observation identification num- 
bers, and the observation date for observations with the 
PACS and SPIRE instruments. These data are available 
from the Herschel Science Archiv^. 

Calibration for the 100 /im and 160 fim PACS bands us- 
ing five standard stars has been found to be good to within 
~ 3% at 100 /zm and 5% at 160 /itm (see PACS Observers' 
Manual). The extended source calibration is good to within 
~ 10% at 70 ^m and 100 ^m and ~ 20% at 160 /xm; we 
adopt the latter numbers here. For SPIRE, the calibration 
is good to within 10% for all bands (see SPIRE Observers' 
Manual). The SPIRE data are in Jy per beam and therefore 
we must know the beam size to convert to a fiux density in 
Jy. Because of diffraction effects, the SPIRE beams are not 
entirely Gaussian - for fiux calculations we use the modi- 
fied beam areas of 423"^ 751"^ and 1587"^ for the 25 ^m, 
350 fim, and 500 ^m bands fsee lSibthoroe et al1l201l[ ). The 
beam area at 250 is ~ 10% larger than would be expected 
for a purely Gaussian beam while those at 350 fim and 
500 fim are approximately that expected. 

We reduce all Herschel data using slightly modified ver- 
sions of the default PACS and SPIRE pipelines within 
the Herschel Interactive Processing Environment (HIPE) 
software, version 7.1. The level 2 pipeline- reduced data 
products from both instruments suffer from striping ar- 
tifacts in the in-scan directions. Additionally, the level 2 
PACS data have artifacts (fiux decrements) around bright 
zones of emission caused by the median filtering baseline 
removal. To limit both the strip ing and the a rtifacts, we 
use the Scanamorphos software] (|Rousselll2012f) . version 9. 
Scanamorphos estimates the true measured value at each 
sky position by exploiting the redundancy in the data. We 
use Scanamorphos without the "Galactic" option; in this 
configuration, the large-scale gradients of the sky emission 
are removed. During a typical Herschel observation, each 
sky position is observed by multiple bolometers in multiple 



scans. These multiple observations of the same sky posi- 
tions can be used to estimate the true sky level. We have 
found Scanamorphos to be very effective at limiting strip- 
ing without the addition of artifacts around bright emission 
zones. Scanamorphos also estimates the photometric uncer- 
tainty at each sky position using the weighted variance of all 
the samples contributing to each pixel. These uncertainty 
maps are useful because at present there is no propagation 
of the errors associated with each processing step in HIPE. 
We use these uncertainty maps in the aperture photometry 
described in [JD 

The astrometry of the PACS and SPIRE Herschel data 
show offsets relative to each other, but also relative to 
higher-resolution Spitzer data. We found that the two 
PACS bands are well-aligned with respect to each other, 
as are the three SPIRE bands. To mitigate the astrometry 
effects, we align the PACS 100 ^m data to the MIPSGAL 
data using point sources detected at both 100 /im and 
24 /xm in the Spitzer MIPSGAL survey (see below). We use 
the same astrometry solution for the PACS 160 /xm data. 
Using point sources detected at 250 /im and 160 /im, we 
then align the SPIRE 250 /im data to the PACS 160 /im 
data and accept the same solution for the SPIRE 350 /im 
and SPIRE 500 /xm data. Offsets were < 10" for all regions, 
in all PACS and SPIRE bands. 

The Herschel data for RCW 120 and Sh 104 w e re firs t 
shown in lAnderson et all (|2010[ ) and iRodon et all (|2010[ ). 
respectively. The data presented in these articles were pro- 
cessed with slightly modified versions of the default PACS 
and SPIRE pipelines within HIPE version 2. . Som e of 
our results d i ffer fr om those of lAnderson et aH (|2010[ ) and 
iRodon et al.l ()2010[ ). in part as a result of this new process- 
ing (see Section IF?^ . 

The new data processing differs from the old processing 
in that: 1) baselines (i.e. drifts on timescales larger than the 
scan leg crossing time) are now estimated by Scanamorphos 
using the redundancy in the data instead of a high-pass 
filter; 2) small-timescale drifts are also derived from the 
redundancy in Scanamorphos (they were either neglected 
[for SPIRE], or attenuated by Fourier filtering [for PACS] 
within HIPE); 3) data projection is done using a projection 
matrix in Scanamorphos for SPIRE data, which prevents bi- 
ases between beam center and the projection center (HIPE 
used a nearest- neighbor projection); and 4) Scanamorphos 
includes relative gain corrections for SPIRE, taking into 
account beam variations from bolometer to bolometer, and 
reducing errors in maps of extended emission; 5) the newer 
version of HIPE used here has updated calibration values. 



3.2. Herschel Hi-Gal 



^ The instrument parameters and calibration for PACS and 
SPIRE listed here are given in the PACS Observers' Manual, 
HERSCHEL-HSC-DOC-0832 Version 2.1 and the SPIRE 
Observers Manual, HERSCHEL-HSC-DOC-0789 Version 2.1. 
Up-to-date versions of both manuals can be found here: 
|http://herschel.esac.esa.int/Documentation.shtml . 
"^ http://herschel.esac.esa.int/Science_Archive.shtml 
^ http://www2.iap.fr/users/roussel/herschel/index.html 



The Hi-Gal survey (jMolinari et al.l 120101 ) . when complete, 
will map the entire Galactic plane using the PACS and 
SPIRE detectors of Herschel. While the Hi-Gal data are 
less sensitive than the observations shown here, Hi-Gal uses 
the PACS detector to observe at 70 /zm, a band not utilized 
with our observation mode. Here we use the Hi-Gal 70 /im 
data for RCW 79, RCW 82, G332. 5-0.1, and RCW 120, the 
four regions that overlap with the current Hi-Gal coverage. 
These data have a spatial FWHM resolution of 6". 
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PACS 








SPIRE 




Name 


Size 


Time 


ObsIDs 


Date 


Size 


Time 


ObsID 


Date 




arcniin. 


sec. 




yyyy-min-dd 


arcmm. 


sec. 




yyyy-mm-dd 


Shl04 


20 X 20 


1730 


1342185575, 1342185576 


2010-10-10 


12 X 12 


629 


1342185535 


2009-10-06 


W5-E 


70 X 70 


13194 


1342191006, 1342191007 


2010-02-23 


65 X 60 


4309 


1342192088 


2010-03-11 


Sh241 


20 X 20 


1730 


1342218725, 1342218726 


2011-04-17 


15 X 15 


770 


1342192093 


2010-03-11 


RCW71 


9x9 


402 


1342189389, 1342189388 


2010-01-16 


6x6 


332 


1342203561 


2010-08-23 


RCW79 


30 X 30 


2768 


1342188880, 1342188881 


2010-01-03 


19 X 19 


837 


1342192054 


2010-03-10 


RCW82 


20 X 20 


1354 


1342188882, 1342188883 


2010-01-03 


12 X 12 


571 


1342192053 


2010-03-10 


G332.5-0.1 


22 X 22 


817 


1342204085, 1342204086 


2010-09-05 


10 X 10 


546 


1342192055 


2010-03-10 


ROW 120 


30 X 30 


3302 


1342185553, 1342185554 


2009-10-09 


22 X 22 


1219 


1342183678 


2009-09-12 



3.3. Spitzer GLIMPSE 

We also use 8.0 /zm data from the Spitzer Galactic Legacy 
Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; 
iBeniamin e^aIll2003^ ■ GLIMPSE e xtends from -f 60° > I > 
-60°, |&! < 1°. The 8.0 ^m IRAC (iFazio et al.ll200l filter 
contains, in addition to the continuum emission of hot dust, 
polycyclic aromatic hydrocarbon (PAH) bands at 7.7 fim 
and 8.6 /im. These molecules emit strongly when excited 
by far-UV photons and thus can be used to trace ioniza- 
tion fronts. The same 8.0 /im IRAC band also shows many 
absorption features, IRDCs, which have been show n to be 
dense molecular clouds (e.g. Simpn et al . 2006; Pillai et al.l 
I2006D . The resolution of the GLIMPSE 8.0 urn data is ~ 2". 



3.4. Spitzer MIPSGAL 



The Spitzer MIPSGAL survey (|Carev et al.ll2009D at 24 /xm 
and 70 fim mapped the Galactic plane from -|-60° < I < 
-60° l&l < 1.0° with the MIPS instrument (i Rieke et al] 
l2004f) . Here we use only data at 24 /im. At 24 /im, the emis- 
sion from Hii regions has two components. First, there is 
emission from very small grains (VSGs) out of thermal equi- 
librium. This e mission is detected in the interior area of th e 
Hii region fsee lWatson et al.ll2008HDeharveng et al.P 2010V 
Secondly, there is thermal emission from the PDR from 
grains that appear to be in thermal equilibrium. These two 
components have rough ly equal fluxes for Galactic bubbles 
(jDeharveng et al.ll2010l ^0| . The resolution of the 24 /im 
MIPSGAL data is - 6". 



3.5. APEX-LABOCA ATLASGAL 

Finally, we use data from the APEX Telescope Large 
Area Survey of th e Galaxy at 870 /im (ATLASGAL; 
ISchuUer et aL|[2009[) . ATLASGAL extends from -^60° > 
I > -60°, |6| < 1.5° and has a spatial FWHM resolution 
of 19". The ATLASGAL survey was performed with the 
Large Apex BOlo nieter CAmera (LA BOCA), a 295-pixel 
bolometer array (jSiringoet al.ll2009D . For RCW120, in- 
stead of ATLASGAL, we use more sensitiv e pointed APEX- 
LABO CA observations first shown in iDeharveng et al.l 
(|2009l) . In the processing of APEX-LABOCA data, large- 
scale low amplitude structures > 2.5' are filtered out. 
The caUbration unc ertainty in ATLASGAL is < 15% 
(|Schuller et al.ll2009D . 

4. The Dust Properties of Galactic Bubbles 

Assuming optically thin emission, the flux from a popula- 
tion of dust grains can be modeled as a modified blackbody: 



where F^, is the fiux density per beam, k^, is the dust opacity 
in cm^ g~^, B^{Td) is the Planck function for dust temper- 
ature Td at frequency v (here in Jysr~^), and A^dust is the 
dust column density in cm^^. The dust opacity law we as- 
sume for the present work is given in Equation [TJ To derive 
the total hydrogen column density and mass (gas-t-dust), 
we must assume a gas to dust ratio, R. The column density 
of gas and dust is: 



iVn-i? 



F„ 



2.8maK^B^iTd)n 



(3) 



where mn is the mass of a hydrogen atom in g, the factor 
of 2.8 accounts for elements heavier than hydrogen, and fi 
is the beam size in steradians. The total mass of dust and 
gas in g m ay then be estima ted from the integrated flux in 
Jy, S^ rcf. lHildebrar^ll983D : 



M = R- 



S,D^ 



K^B„{Td) 



(4) 



F^ OC K^ B^{Td) Nd^st : 



(2) 



where D is the distance in cm. The preceeding equations as- 
sume that the measured flux at frequency v is entirely due 
to thermal dust emission. For the present work, at 350 /im, 
w e assume the functiona l form for the dust opacity given 
in lBeckwith et al.l (|1990D . k„ = lO(i^)'', where z/ is in THz 
and we assume /3 = 2. This leads to an opacity at 350 /im of 
^350 = "7.3 cm^ g^^, which we use throughout. For all mass 
calculations, we assume that the dust-to-gas mass ratio has 
a value of 100. These values for k^, and R are uncertain, 
which leads to large uncertaint ies in the mass a nd co lumn 
density estimates. For example. iPreibisch et al.l (|1993[ ) esti- 
mate that using the above formulation to compute masses 
results in uncertainties of up to a factor of five. Errors in 
distances for some sources may be up to a factor of two 
in extreme cases, further increasing the uncertainty of the 
mass estimates. 



4.1. Aperture Photometry 

Aperture photometry allows us to derive the average dust 
properties within an aperture of arbitrary size and shape. 
The measured flux within an aperture contains the com- 
bined emission along the line of sight. Most of the regions 
studied here are optically visible and thus the extinction 
along the line of sight is relatively low. The contribution to 
the fluxes from dust along the line of sight not associated 
with the H ii region is therefore minimized, but cannot be 
completely absent. 

We deflne the apertures by eye to include contiguous 
regions of dust emission from a volume that can be charac- 
terized by a single temperature value. This process is sub- 
jective; we use the three-color images shown in Figure [1] to 
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Fig. 2: Same as Figure [I] but showing the regions used in the aperture photometry. As in Figure [TJ the three-color images 
are composed of Herschel data for each source with 500 fim in red, 250 /im in green and 100 fiui in blue. The apertures 
used for the aperture photometry measurements are shown in gray outline. 



help define regions where the temperature variation within 
the region is not severe. In addition to specific regions of 
interest, we define for each Hii region a large aperture that 
includes all emission associated with the Hii region itself, 
but not the emission associated with local clouds or fila- 
ments. The apertures are shown in Figure [2] on top of the 
same Herschel data shown in Figure [T] In total, we define 
129 apertures. 

We visually categorize the apertures into five classi- 
fications based on their location and appearance in the 
Herschel data. The five aperture classes are: "Entire" for 
the largest apertures that include all the emission from the 
Hii region (including that of hte FDR), "Other Hii" for 
the apertures that include emission from other H ii regions 
in the fields, "PDR" for locations along the PDRs, "Point 



source" for apertures containing only the emission from 
sources that appear point-like at 100 /im (including pos- 
sibly their outer envelope detected at longer wavelengths), 
and "Filament" for apertures that enclose extended struc- 
tures detected at 250 /xm that are not along PDRs and that 
are not the longer- wavelength emission of point-like objects 
detected at 100 /im. Because filaments seen in absorption at 
8.0 /im must be dense and therefore may be especially cold, 
we subdivide the "Filament" class into those seen in ab- 
sorption at 8.0 /xm (IRDCs), and those that are not seen in 
absorption at 8.0 /im. Excluding the "Entire" and "Other 
Hii ' classes, which both contain eight apertures, there are 
113 apertures. Of these, 49 are classified as "PDR", 46 as 
"Filament" (with 17 as infrared dark filaments), and 18 as 
"Point source" . 
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Fig. 2: - continued. 



Many of the ap ertures for RCW 120 an d Sh -104 were de- 
fincd previously in I Anderson et aLl (J201 0') and Rodo n et al.l 
(.2010.) . respectively. We have made sonic modifications to 
the apertures shown in these works: wc have changed the 
shape of some apertures, added new apertures to highhght 
additional regions of interest, removed apertures that ei- 
ther cannot be well-characterized by a single dust tem- 
perature, and given some apertures a different name. The 
new aperture shapes should better contain dust that can be 
characterized by a single temperature, and are more consis- 
tant between all regions in the current sample. We removed 
apertures in the direction of the "interior" of the bubbles 
because of the difficulty of fitting a single dust temperature. 
The aperture names used here better reflect the location of 
the aperture and again are more consistant with the naming 
convention employed for the other regions in our sample. 



In addition to the "bubble" H ii regions that we targeted 
in our observations, there are eight additional sources in 
the Herschel fields that have characteristics of Hii regions. 
All of these eight sources have spatially coincident IR (from 
Herschel and Spitzer) and radio emission from the NVSS or 
SUMSS. Together spatially coin cident IR and radio emis- 
sion point to a therrnal source ([Haslam fcOsbornd 119871 : 
iBroadbent et al.iri989HAnderson et al.ll201l[ ). e.g., a plane- 
tary nebula (PN) or an H ii region. 

lAnderson et al.l (J2012t ) derive IR color criteria for dis- 
tinguishing between PNe and Hii regions. To determine 
the classification of the eight other Hii regions, in addi- 
tion to the aperture photometry described below we com- 
pute aperture photometry using the Wide-field Infrared 
Survey Explorer (WISE) data at 12 /xm and 22 /im. All 
the Hii regions tested satisfy the criterion found by 
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lAnderson efaP (|2012D : Iog(i^i6o/i^i2) > 1-3, where Fieo 
is the flux at 160 fiui from Herschel and F12 is the flux at 
12 /im from WI SE and in fact all have log(Fi6o/-Fi2) > 2. 
lAnderson et al.l (|2012f ) found that 10% of their sample 
of PNe satisfy Iog(i^i6o/-Fi2) > 1-3 and no PNe satisfy 
log(Fi6o/-Fi2) > 2. Together with the spatially coincident 
IR and radio emission, these criteria suggest that these 
sources are bona fide Hii regions. There are no WISE data 
available for the additional H 11 region located on the PDR 
of Sh2-104. Since this source satisfies the criteria for UC H 11 
regions in IWood & ChurchwcU (1989a]) and has detected 
CS emission ([Bronfman et al.lll996l ). it is almost surely an 
H II region and not a PN. We assume throughout following 
that the eight other H 11 regions are at the same distance as 
the bubble H 11 regions that were targeted. This assumption 
may not hold for H 11 regions well-separated in angle from 
the bubble Hii regions, i.e., for one of the Hii regions in 
the field of G332 and the second H 11 region in the field of 
Sh241. 

Because our data have a range of resolutions, we smooth 
all data to the lowest resolution SPIRE 500 /im data, which 
has a resolution of 37". To mediate pixel edge effects, we 
rebin all data sets so they have the same pixel size and 
location using the Montage softwareQ. 

For each H 11 region, we compute the fiux at all available 
wavelengths for each aperture. In this process we utilize the 
Kang softwarqj, version 1.3, which contains routines to per- 
form aperture photometry with apertures of arbitrary size 
and shape, and accounts for fractional pixels. We estimate 
the errors at each wavelength by summing in quadrature 
the calibration uncertainties given in SJ3] and the photomet- 
ric error as calculated by Scanamorphos. For Herschel data, 
we compute the photometric error by adding in quadrature 
the individual errors of all pixels within the aperture from 
the error maps produced by Scanamorphos. For the other 
data sets, we estimate the photometric error using the stan- 
dard deviation of the background region multiplied by the 
square root of the number of pixels in the aperture. 

We subtract from each aperture fiux, at each wave- 
length, a background value. For apertures at the center of 
the frame (mainly the "Entire" , "Other H 11," and "PDR" 
classes) we use a single background value at each wave- 
length for all apertures. Because apertures in these clas- 
sifications generally have high fiuxes, we have found that 
the exact choice of background value does not severely im- 
pact our results. A single background value, however, does 
not approximate the true background level accross the en- 
tire field. Therefore, for apertures located away from the 
center of the frame we define a local background value. 
The removal of a background is necessary because the true 
"zero-point" for our observations is unknown. Perhaps more 
importantly, however, there are potentially contributions 
from multiple emission regions along the line of sight. In 
the Scanamorphos data reduction, we remove the general 
gradients of the sky emission within the field of view (by not 
using the "Galactic" option) and therefore we believe the 
fiat background approximation is valid; attempts to remove 
a more complicated background were met with difficulty 
due to the relatively small size of our fields. 

We construct a spectral energy distribution for each 
aperture using the background-corrected fiuxes at wave- 
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Fig. 3: Example SEDs fits for regions within W5-E (top) 
and RCW79 (bottom). The solid curves show the /? = 2.0 
trial and the dashed curves show the /3-free trial. We do 
not have 70 fim or 870 /im data for W5-E, and we exclude 
the 870 /im point from the fit to the "Entire" aperture for 
RCW79. 



lengths from 70 /im to 870 /im. Using the "MPFIT"0 least- 
squares minimization routines in IDL (Markwardt 20091 ), 
we then fit Equation ^ in two trials: once with the nor- 
malization and dust temperature as free parameters and /? 
fixed to a value of 2.0 and once with the normalization, dust 
temperature, and /3 as free parameters. In the SED fits, we 
generally include data at both 70 /im and 100 /im, although 
there is likely some contribution from a warmer dust com- 
ponent at these wavelengths. We do, however, exclude data 
at 70 /im from IRDC apertures because at this wavelength 
they may be optically thick and we exclude 870 /im data 
from the "Entire" apertures because these data are not sen- 
sitive to large-scale diffuse emission. Emission from warmer 
dust is detected mainly toward the interior of the bubbles 
in our sample - it has minimal impact on the parameters 
derived here (see tj5.4p . Using Equation 01 we also calculate 
the mass within each aperture. 

The measured flux within a photometric bandpass is a 
function of the shape of the SED across the bandpass and 
of the shape of the bandpass itself. Each measured fiux is 
therefore a function of temperature and /3. To correct for 
this issue, we perform a "color correction" on the measured 
fiuxes. During the spectral fitting, we iteratively apply color 
correction factors by repeatedly fitting a temperature and 
/3 then applying the corresponding correction factors until 



see |http://purl. com/net /mpfit| 
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successive fit results are unchanged. For SPIRE, these color 
correction factors, given in the SPIRE Observers' Manual, 
are < 10% for all values of /3 conside red here. For PACS, th e 
color correction factors are given in iPoglitsch et al.l (J2010l) . 
The PACS color correction factors can be quite large at low 
temperatures, especially for the 70 ^m band. For example, 
at 10 K, the correction factor is 3.65 for the 70 /im band and 
1.71 for the 100 /im band. At 15 K, the respective factors 
are 1.61 and 1.16. While these correction factors are large, 
we find that the net impact on the derived temperature is 
small. For the apertures used here, the mean absolute tem- 
perature difference between apertures fit without a color 
correction and the same apertures fit with a color correc- 
tion is 0.6 K, or -3%. 

The results of the aperture photometry analysis are 
given in Table 3. This table lists for each aperture the name 
shown in Figure [H mean Galactic longitude and latitude, 
the angular size of the aperture in square arcminutes, the 
derived dust temperature and mass for the two fit trials, 
the /3-values for the /3-free trial, and the classification. The 
errors in dust temperature and /3 account for calibration 
uncertainty and photometric errors. We show in ii5.2l that 
they are good estimates of the errors in the derived quanti- 
ties. Example fits are shown in Figure [3] for three apertures 
within W5-E and RCW79. 



4.2. Maps of Dust Properties 

The high-angular resolution Herschel data also allow us to 
produce two-dimensional maps showing the distributions 
of temperature and column-density. With such maps, we 
may examine both small- and large-scale variations in dust 
properties. These maps have the resolution of the 500 /xm 
data, 37" and are available from the authors upon request. 



4.2.1. Temperature Maps 

We construct temperature maps for each of the eight re- 
gions in our sample by fitting the SED extracted at a grid 
of locations using Equation [2] When doing so, we use the 
same chi-squared minimization routine used for the aper- 
ture photometry and again leave the column density as a 
free parameter. In the creation of the temperature maps, we 
again include data at 70 /im; including the 70 /x m emission 
differ s from the method of other authors (e.g., iHill et aLl 
1201 ih . The effect of including data at this wavelength is 
explored in H5.4I 

While the wavelength coverage of Herschel and the num- 
ber of data points is theoretically sufficient to simultane- 
ously fit for Td and /3, we find that the uncertainties in 
such fits are large and the correlation between the derived 
temperatures and the values derived in the aperture pho- 
tometry is weak. We therefore hold (3 fixed to a value of 2.0 
when creating the maps, which is the average value found 
in the aperture photometry. 

We remove a flat average background level from the 
smoothed, rebinned data at each wavelength. The back- 
ground value was defined such that 99.99% of the pixels in 
the resultant images have an intensity greater than zero. 
This ensures that nearly all locations have positive fiux 
values. This method of background subtraction is different 
from the method of the aperture photometry but is nec- 
essary to ensure that background-subtracted pixels have a 



positive intensity value, and thus that we will be able to 
determine the dust temperature and column density over 
the entire map. 

We show temperature maps for the eight regions in 
Figure E] for the same map areas shown in Figures [1] and 
[21 The highest temperatures, ~ 40 K, are found in the di- 
rection of the interior regions of the bubbles and toward 
other local Hii regions. The coldest temperatures, '^ 12K, 
are found for the IRDCs and other regions designated as 
"filaments" in the aperture photometry (t j4.ip . Errors in 
the derived temperatures for individual pixels are generally 
less than 10%. Large patches of very high temperatures on 
the periphery are in general not real, but are caused by low 
emission (e.g., to the north-east of Shl04, to the north of 
Sh241, and to the north-east of RCW82). 

Figure [4] shows a wide range of temperature struc- 
tures. For Sh 104, Sh 241, RCW 79, RCW 120, and RCW 82, 
the bubble interiors show increased temperatures. These 
temperatures are quite uncertain though and individual 
SEDs often show contributions from multiple components 
(see t )5.4|l . The hottest temperatures for RCW 71, however, 
are in the PDR. For G332.5-0.1 the hottest dust tem- 
perature is also along the PDR, to the north, coincident 
with the UC and compact Hii regions. W5-E shows lit- 
tle variation in temperature although there are warmer lo- 
cations associated with the other Hii regions in the field. 
Variations within the fields of individual regions are di- 
cussed in Appendix \K\ 

The temperature map values agree with the results of 
the aperture photometry, but there are large differences 
for colder regions. To compare the two methods of deriv- 
ing tempatures, we compute the average dust temperature 
found in the temperature map within the apertures used 
previously. In Figure \5\ we plot the temperature derived 
in the aperture photometry versus the average tempera- 
ture map temperature for each aperture. We find that the 
two temperatures are on average 12% different, or roughly 
2.3 K. This difference is pronounced for apertures with tem- 
peratures less than 15 K and warmer than 30 K, for which 
the temperature map values are sytematically higher or 
lower, respectively. These discrepancies are due to the dif- 
ferences in background treatment between the two meth- 
ods. Because we evaluate the background individually for 
each aperture, we believe the aperture photometry dust 
temperature values are more accurate. The circled points in 
Figure [5] are IRDCs in the field of RCW 120; they account 
for much of the difference at low temperature. The differ- 
ence in temperature for these cold regions is due to our 
assumption of a fiat background in the temperature maps. 
For cold regions this assumption causes us to overestimate 
the temperature, while the inverse is true for warm regions. 
We do not expect a perfect correlation because the aperture 
photometry temperatures are weighted toward the highest 
emitting regions within the aperture. The mean standard 
deviation of temperature map values within the apertures 
is 2. OK. 



4.2.2. Column Density Maps 

From the temperature maps, we also create column density 
maps using Equation [31 As mentioned above, we evaluate 
Ki, at 350 fim, and assume a value of K350 — 7.3cm^g~^ 
and a dust-to-gas mass ratio of 100. At each grid loca- 
tion in the column density map, we calculate the Planck 
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Fig. 4: Temperature maps created by fitting the SEDs extracted at a grid of locations. All panels have the same color 
scale and range from 15 K to 35 K. The green curves show the approximate extent of the bubble interiors. 



function value By{Td) using the derived temperature from 
the temperature map. For the flux at 350 /Ltm, we use the 
Herschel SPIRE image data, after removing a background 
value defined as before such that 99.99% of the pixels at 
350 /zm have an intensity greater than zero. The removal 
of a background value produces a more realistic picture of 
the column density associated with each region, rather than 
the total integrated column density along the line of sight. 
The resultant maps are shown in Figure [5] 

The column density distributions in Figure [B] shows en- 
hancements along the bubble PDRs. These enhancements 
likely are from the material swept up during the expansion 
of the Hii region. For RCW71, however, the PDR is barely 
enhanced relative to the background. Other prominent col- 
umn density enhancements in the fields are generally caused 
by filaments, especially IRDC s, and unresolved sourc es. For 
comparison, the IRDCs in .Peretto fc Fulleii (|2010( ) have 
typical column densities of > lO^^cm"^. The highest col- 
umn densities are generally associated with condensations 
with detected protostars inside (i.e. Condensation 1 in 
RCW120). Column density values toward the bubble in- 
teriors are uncertain due to uncertainties in temperature 



and because of the low flux values found there. The spu- 
rious temperature values near the image edges mentioned 
previously result in spurious column density values as well 
(e.g., south-west of Sh241 and south-east of RCW71). A 
more complete discussion for individual regions is given in 
the Appendix. 

PDRs and cold filaments share nearly the same mean 
column density values. We calculate the mean column den- 
sity value shown in Figure [6] within the previously defined 
apertures for regions in the "PDR" and "Filament" classes. 
The average value for the "PDR" class is 1.8 x lO^^cm"^ 
while it is 1.9 x 10^-^ cm~^ for the "Filament" class. The 
above average values were calculated after first taking the 
log of the column density values, so large values of column 
density would not bias the result. 



5. Discussion 

5.1. Comparison of Aperture Classes 

With the calculated values of /3 and Td in Table 3, we may 
investigate how the dust properties change with location in 
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Fig. 4: - continued. 



the Hii region environment using the aperture classifica- 
tions. The easiest way to see differences in dust properties 
is in the shapes of the SEDs themselves. We show average 
SEDs for the classifications in Figure [71 To create Figure [71 
we normalized the SED for each aperture by the highest flux 
and averaged all apertures for a given classification. The 
curves are therefore unweighted averages of the normalized 
fiuxes. The "Entire," "Other Hii," and "PDR," classes have 
similar shapes, as do the "Filament" and "IRDC" classes. 
The "PS" class has a mean SED shape in between the two 
groupings. The similar shapes of the SEDs hint at similar 
dust properties for these classifications. Compared to the 
other classifications, the "Filament" and "IRDC" classes 
peak at longer wavelengths and show more emission in the 
FIR part of the SED. This shows that these apertures have 
colder mean temperatures. 

We list the mean dust properties for each aperture clas- 
sification in Table 4, which shows the average dust tem- 
perature and /^-values for the /3 — 2 trial and the /3-free 
trial, together with their l-cr standard deviations, the per- 
centage of apertures that have peak fluxes long-ward of 



100 /im, the percentage of apertures that have peak fluxes 




15 20 25 30 35 

Aperture Temperature (K) 

Fig. 5: The temperature derived using aperture photom- 
etry versus the average temperature in the temperature 
map derived using the same apertures. The line is not a 
flt but rather shows a 1:1 relationship. The circled sources 
are IRDCs from the held of RCW120; they account for 
much of the difference at low temperature. 
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Fig. 6: Column density maps created using Herschel 350 /im data and the temperatures shown in Figure |4l The scaling 
for all panels is the same; the values range from lO^*' cm~^ to lO^^*^ cm^^. The green curves show the approximate extent 
of the bubble interiors. Enhancements are seen along the bubble PDRs. 



long-ward of 160 /im, and the number of apertures of a 
given classification. The seventh and eigth columns deal- 
ing with the wavelength of the peak flux provide a check 
against the results of the temperature fitting - SEDs that 
peak long-ward of 160 /im have dust temperatures of < 12 K 
while those that peak long- ward of 100 /im have dust tem- 
peratures of < 21 K, assuming /3 — 2. The row labeled 
"All Filaments" includes filaments of aperture classifica- 
tion "F" and infrared dark filaments of aperture classifi- 
cation "F (ID)" from Table 3. The rows labelled "IRDC" 
and "Not IRDC" contain only data from inner-Galaxy H ii 



regions for which we have 8.0 nm data (RCW79, RCW82, 
G332. 5-0.1, and RCW120). 



in temperature has implications for 
formation. The Jeans mass scales as 



A difference 
subsequent star 
Mj oc T'^/^ /9~^/^, where p is the mass volume density, and 
therefore higher temperatures lead to higher values of Mj. 
The average temperature of the PDR classification, 25.6 K, 
is significantly higher than that of the filaments, 16.9 K. 
For the above average temperatures, if the gas and dust 
share the same temperature and the density in the PDR 
and the filaments is comparable, the Jeans mass would be 
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Fig. 6: - continued. 



^ 80% higher in the H ii region PDRs compared to the fila- 
ments. This would imply that on average higher mass stars 
will form in the H ii region PDRs through gravitational col- 
lapse compared to the filaments. We stress, however, that 
the density is a large unknown here and further observa- 
tions are needed to better constrain this parameter. 

Table 3 and Figure [7] show that there are two broad clas- 
sifications of apertures that can be grouped by their similar 
temperature values. The first is a "warm" group character- 
ized by temperatures of ~ 26 K that includes the "Entire" , 
"Other Hii", and "FDR" aperture classifications. There 
is no statistical difference between all possible combina- 



tions of aperture classifications within the "warm" group, 
as shown by a Kolmogorov-Smirnov (K-S) test which as- 
sesses the probability that two distributions belong to the 
same parent distribution. The second "cold" group consists 
of apertures classified as "Filament" or "IRDC" and is char- 
acterized by dust temperatures of ~ 17 K. Some apertures 
in the "PS" classification have dust properties similar to 
those in the "warm" group, while some have dust proper- 
ties more similar to those in the "cold" group; we do not in- 
clude "PS" apertures in either group. This likely refiects the 
range of evolutionary states in this classification: the colder 
apertures of the "PS" class are in an earlier evolutionary 
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Fig. 7: The average normalized SEDs for the aperture classes. Results from the "PDR" class are shown in dark gray while 
those in the "Filament" class are shown in light gray. Where the two "PDR" and "Filament" distributions overlap, we 
show light gray with dark gray hatching. Lines show the average normalized SEDs for the "Entire" , Other "Hii," "PS," 
and "IRDC" classes; the "IRDC" class SED is similar to that of the "Filament" class while the SEDs of the "Entire" , 
and "Other Hii" classes are more similar to that of the "PDR" class. The "PS" class has a mean SED shape in between 
that of the "PDR" and "Filament" classes. The points at 8.0 /im and 24 /im are shown here although these data were 
not used in the fits. 



state than the warmer ones that may be heated by a cen- 
tral young stellar object. In Figure [51 we show graphically 
how the "warm" and "cold" aperture groups, which were 
defined according to their aperture classifications and not 
according to their dust temperatures, have different tem- 
perature distributions. Roughly half of the apertures in the 
"warm" group have an emission peak long-ward of 100 /xm 
while all the apertures in the "cold" group meet this cri- 
terion. There are no apertures in the "warm" group whose 
wavelength of peak fiux emission is long- ward of 160 /im, 
while 13% of the apertures in the "cold" group meet this 
criterion. 

Filaments that are infrared dark at 8.0 /im are colder 
than those that are not infrared dark at 8.0 /<m. When /3 
is held fixed to 2.0, the "IRDC" filaments have a mean 
temperature of 15.5 ± 3.5 K while the "Not IRDC" fila- 
ments have a mean temperature of 18.8±2.9K; these values 
change to 14.7 ± 3.6 K and 17.8 ± 2.2 K respectively when 
/3 is allowed to vary. A K-S test shows that the tempera- 
tures derived for the two groups when /3 was held to 2.0 are 
significantly different. This suggests that IRDCs are colder 
versions of the filaments we identified. The detection of an 
IRDC relies on a strong IR background that the cloud can 
absorb against and some of the filaments not seen in ab- 
sorption may simply be lacking sufficient background. Cold 
dust temperatures Td < 12 K are a lso found in less mas- 
sive filaments of nearby c omplexes (jBontemps et al.ll2010t 
lArzoumanian et all 1201 11) and thus such temperatures are 
not restricted to the more massive Galactic IRDCs. 



As shown in Table 4, apertures in the "warm" and 
"cold" groups share a similar distribution of /3 values. 
We show these distributions graphically in Figure [HI Both 
groups are centered about /3 = 2.0. The average error in /3 
is 0.3, which likely causes some of the spread in the distri- 
bution. 



5.2. The (3 — Td Relationship 

As discussed in the Introduction, there is significant ob- 
servational evidence for an anti-correlation between f5 and 
dust temperature. It is therefore reasonable to expect the 
same correlation for our data. In Figure[TU]we plot the dust 
temperature against (3 for all apertures. We do not find a 
strong relationship between the two quantities. The data 
points in Figure [TU] are coded such that apertures in the 
"warm" group are shown with filled circles and those in the 
"cold" group are shown with open circles. The "PS" aper- 
tures are shown wit h star symbols. The tr e nd lines found by 
iDupac et aP (2003"), 'Desert et all (|2008[ ). IVenezian i et all 
i201GD , and [Paradis et al. (2010|) are over-plotted, as is the 
best-fit line to our data. The trend lines found by previous 
others do not agree with our trend line. Most of these pre- 
vious studies have used different wavelengths to those used 
here. 

Using the same functional form as lDesert et al.l (|2008[ ). 
our best-fit line follows the relationship: 



P = (4.3 ± 0.5) t;^ 



-0.2±0.04 



(5) 
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Table 4: Dust properties by aperture classification 



Class. 



13 = 
<Ta> 
K 



K 



<Td> 
K 



/3 free 
aTd <I3> 
K 



•y/i 



> 100 /im 

% 



> 160 /Jin 

% 



Entire 


25.3 


2.4 


26.5 


2.9 


1.9 


0.3 


50 





8 


Other HII 


26.1 


2.0 


26.2 


3.2 


2.0 


0.3 


37 





8 


PS 


21.3 


3.3 


22.2 


3.1 


1.8 


0.2 


83 





18 


PDR 


25.8 


4.0 


24.4 


3.3 


2.1 


0.4 


61 





49 


All Filaments 


16.8 


3.0 


16.8 
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2.1 


0.5 


100 


13 


46 


IRDC 


15.5 


3.5 


14.7 


3.6 


2.2 


0.4 


100 


29 
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Not IRDC 


18.8 


2.9 


17.8 


2.2 


2.2 


0.2 


100 
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Fig. 8: Dust temperatures of warm and cold aperture groups. Shown are the distributions for the apertures in the cold 
group ( "Filament" and "IRDC" classifications) in light gray and the apertures in the warm group ( "Entire," "Other H ii," 
and "PDR" classifications) in dark gray. We do not include in this figure apertures with classification "PS" . Where the 
two distributions overlap, we show light gray with dark gray hatching. The left panel shows dust temperatures derived 
when /3 was held fixed to a value of 2.0 while the right panel shows the dust temperatures derived when /3 was allowed 
to vary. It is evident that these two aperture groups have different dust temperature distributions. 



A linear fit follows the relationship: 
/3 = 2.6 ± 0.1 - (0.02 ± 0.005) Td 



(6) 



Both fits are rather poor and share the same reduced x^ 
value of 2.2, with the same number of free parameters. The 
linear fit in Equation [B] shows essentially no relation be- 
tween the two quantities. Given that the linear relationship 
is the simpler of the two forms, and both have the same sta- 
tistical significance, we believe it is more representative of 
the data shown here. Amongst the other studies considered. 
Equation E] shows the best agreement with the relation of 
iDupacet al1(|2003D . 

It is interesting that the cold and warm groups alone 
appear more consistent with a /3 — T^ relationship than 
when they are combined. A fit of the form in Equation [5] 
to the cold group has a reduced x^ value of 1.4, while a fit 
to the warm group using the same functional form has a 
reduced x^ value of 2.2. These fits, while better, are still 
rather poor. We defer further discussion of this point to a 
future publication. 

5.2.1. Simulations of the f3 — Td Relationship 

We perform simulations to better understand the effect of 
calibration and photometric uncertainties on the derived 



Td and /3-values. According to their user manuals, the flux 
calibration for the PACS and SPIRE instruments is stable 
in time and the calibration for the photometric bands of a 
given detector is correlated. We model four different combi- 
nations of calibration errors: -f 5% for all PACS bands and 
+5% for all SPIRE bands, -1-5% for PACS and -5% for 
SPIRE, -5% for PACS and -1-5% for SPIRE, and -5% for 
PACS and -5% for SPIRE. We use +5% for ah 870 fiui cal- 
ibration errors. We chose these values for the uncertainty 
estimates to assess how different combinations of uncertain- 
ties may affect the derivation of the f3 — T relationship. The 
values are not necessarily consistant with the instrumental 
and photo metric unc ertainties discussed in more detail be- 
low. Ander son et al.l |2012) found using Herschel data that 
the photometric uncertainties are in most cases dominated 
by the choice of background region. We did not rigorously 
estimate this source of uncertainty here but we assume tha t 
the values are similar to those in [Anderson et al.l (|2012[ ) . 
For the photometric uncertainty, we take the average un- 
certainty they give for their photometry of H ii regions: 7%, 
10%, 11%, 13%, and 15% for the 70 ^m, 160 ^im, 250 /xm, 
350 /xm, and 500 ^m bands. We assume the photometric 
uncertainty at 100 /xm is 10% and the error at 870 /im is 
15%. 

We simulate the emission from dust that has a range of 
temperatures from 10 to 35 K in steps of 5 K, in a manner 
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Fig. 9: Dust opacity spectral index /3 of warm and cold 
aperture groups. Shown are the distributions for the aper- 
tures in the cold group ("Filament" and "IRDC" classifi- 
cations) in light gray and the apertures in the warm group 
("Entire," "Other Hii," and "PDR" classifications) in dark 
gray. Where the two distributions overlap, we show light 
gray with dark gray hatching. The two aperture groups 
have only slightly different dust emissivity index distribu- 
tions. The typical error in /3 for individual apertures is 
-0.3. 
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Fig. 10: The jS-Td relationship derived from the aperture 
photometry. Data points from the "warm" group are shown 
with filled circles while those from the "cold" group are 
shown with open circles. Data points from "PS" aper- 
tures are shown with star symbols. We plot on top of the 
da ta the best fit line a s well as the best fi ts lines found 
bvlDupac et al.l (l2003l). iDesert et a ll (I2008D . IP aradis et al.l 
(|2010D . and lVeneziani et all (|2010t ). 



similar to that of IShettv et al.l (|2009al ). For this simulation 
we use a single /3 value of 2.0. For each (3 — Td combination, 
we calculate the fluxes of a simulated gray-body SED. We 
then randomly extract values from Gaussian photometric 
error distributions with Icr-levels given above. We add the 
calibration uncertianties in quadrature to the photometric 



errors. The resultant synthetic SEDs have flux values var- 
ied by representative calibration and photometric uncer- 
tainties. We have 16 trials from the combinations of four 
wavelength options and four calibration options. 

We show the results of this modeling in Figure [11] for 
the four trials with calibration errors of -1-5% for PACS and 
—5% for SPIRE. Each data point in Figure [TT] represents 
the derived Td and /3-values for a single simulated SED; 
there are 600 data points in total, 100 for each input tem- 
perature. Points sharing a common input temperature also 
share a common shade of gray. All four trials with differ- 
ence calibration errors that share a wavelength combination 
have results similar to that shown in Figure [11] although 
the scatter in the data points is generally larger when the 
PACS and SPIRE calibration errors have inverse signs. 

A number of points are evident in Figure 1111 First, for 
each (3 — Td pair, a false relationship is produced. Secondly, 
a large number of wavelengths helps greatly in reducing the 
spread in the calculated values and the uncertainty in the 
derived flt. Excluding either the 70 /im or the 870 /im data 
results in ~ 50% larger flt uncertainties; excluding both 
wavelengths results in ~ 100% larger fit uncertainties in 
both /? and Td. Finally, the most likely value found is that 
of the input value; there is no systematic offset between the 
derived and input j5 — Td pair. 

Could the f3 — Td relationship be falsely produced by our 
analysis? To better assess the effect of the Herschel calibra- 
tion and photometric uncertainties on our derived Td and 
/3-values, we repeat the simulation above using as input 
temperatures the 129 temperature values from the /3- fixed 
trial derived in the aperture photometry (see Table 3), and 
/3 = 2.0. When constructing the simulated SEDs, we in- 
clude data at the same wavelengths as were included in our 
129 aperture photometry fits: 42 data points with 70 /im 
and 870 /zm data, 2 data points with 70 /im but no 870 /im 
data, 16 data points with 870 /im data but no 70 /im data, 
and 69 data points with neither 70 /im nor 870 /im data. 
We use for the calibration errors one value for all PACS 
bands and one value for all SPIRE bands. This calibration 
error is drawn randomly from a uniform distribution, us- 
ing as a maximum the percentage error given previously in 
Section [31 We treat the photometric errors in the same way 
as described above and again add these two sources of error 
in quadrature. We simulate the 129 SEDs 1,000 times, us- 
ing for each of the 1,000 trials different calibration errors. 
We fit both a linear regression line and also a power law 
regression to each of the 1, 000 sets of 129 points. 

In Figure [T^ we show an example from one of the 1, 000 
trials: the trial whose linear regression fit parameters are 
the closest to the median fit parameters from all 1, 000 tri- 
als. As in Figurc[TTl each data point in Figure [T^ is the fit to 
a simulated SED that includes calibration and photometric 
errors. If there were no sources of error added, we would re- 
cover the 129 derived temperatures, all with /3 = 2.0. The 
simulated temperature and /3- values shown in Figure IT^ are 
visually similar to our derived values shown in Figure 1101 
The median linear fit is: 



;9 = 2.6 ± 0.2 - (0.02 ± 0.01) Td , 
and the median power-law fit is 

13= (4.o±i.3)r7°-2±°-^ 



(7) 



(8) 
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The simulated fits again show a weak inverse relationship 
between /? and T^, with the same values (within the errors) 
as what was found in Equation [S] and Equation [51 

From Figures [TU] and [1^ there appears to be less varia- 
tion in /3 in the simulated data than in the real data. The 
summed square of the difference in /3 between the linear fit 
and values of /3 in Figure [TU] is 15.6. The median summed 
squared difference in the 1, 000 fits to the synthetic data de- 
scribed above is 13.0 with a large standard deviation of 8.2. 
We therefore conclude that most of the variation in /3 can 
be explained by calibration and photometric uncertainty 
but that there may also be variations in /3 not explained by 
these sources of uncertainty. 

5.2.2. Comparison with Previous Results Using These Data 

The lack of a strong relations hip between Td and is in 
contrast to the su ggestion of [Anderson et al.l (|2010f ) and 
iRodon et al.l (|2010f ) . Using an earlier processing of the same 
data employed here for RCW120 and Sh-104 (except for 
the Herschel Hi-Gal 70 fj.m data for RCW120), these au- 
thors' results were consistent with the existence of a rela- 
tionship between Td and /3— T, although these works did not 
consider the error simulations that we perform here. The 
methods between their works and ours are nearly identical 
and therefore it seems likely that the updated processing 
of the Herschel data has caused some of the differences. 
The median absolute flux differ ences between o u r proc ess- 
ing and the processing used in [Anderson et al.l ([20101 ) for 
RCW 120 (after rebinning to the same pixel grid and sub- 
tracting a background such that 99.99% of the pixels have 
positive values) are 2.3%, 4.7%, 0.7%, 0.7%, and 0.8% for 
the 100 /xm, 160 /im, 250 /im, 350 fim, and 500 fim data. 
The median absolute flux differences for Sh-104 are similar: 
0.8%, 1.2%, 0.4%, 0.3%, and 0.5% for the 100 fim, 160 /.tm, 
250 /zm, 350 /im, and 500 fxui data|3 

The reprocessing has affected some derived values of 
Tdust and P, but these differences are generally within 
the previously stated uncertainties. Because w e use Hi-Gal 
TOum data here in plac e of the MIPSGAL ([Carev et al.l 
[2009D 70 /im data used in [Anderson et al.[ ([2010t ) the differ- 
ences in derived Td and /3- values reflect the combined effect 
of the new processing (including changes to the calibra- 
tion) and the different 70 /im data sets. For the trial when 
/3 was allowed to vary, the median absolute temperature 
difference is 2.4 K. About half of the apertures in common 
have temperatu r e valu es within the uncertainties given in 
[Anderson et al.[ ([2010( ). The median absolute difference in 
/3 is 0.3, and over h alf of these differences a re within the un- 
certainties given in [Anderson et al.1 ([2010[). The da t a use d 
here for Sh-104 are the same as that of [Rod6n et al.l ([20101) . 
albiet reprocessed. We find that, for the trial when /? was 
free to vary, the median absolute difference in temperature 
betwe en what we find here and that shown in [Rod6n et al.l 
([2010t) is 1.2 K; in all cases this median absolute tempera- 



^ We also checked the flux differences between the current 
Herschel processing data and a processing employing HIPE ver- 
sion 4.0 and Scanamorphos version 2.0 and find very similar re- 
sults for both regions. Furthermore, the differences between the 
data processed with HIPE version 4.0 and the data used here 
are < 0.3% for all photometric bands for both regions. This indi- 
cates that the current processing is stable and we do not expect 
future versions of HIPE or Scanamorphos to signiflcantly alter 
the flux values. 
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Fig. 11: Fits of Td and (3 using simulated data. Each data 
point is the result of an SED fit to simulated data account- 
ing for calibration and photometric uncertainties. The sim- 
ulated SEDs were created with input temperature values 
from 10 to 35 K in steps of 5 K and (3 = 2.0. Points with a 
common input temperature share a common shade of gray. 
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Fig. 12: Fits of Td and /? using simulated data, similar to 
Figure [TlJ We used the fitted temperatures from Table 3 
to create the simulated data, and /3 — 2.0. Shown is one 
of 1,000 trials using different simulated calibration errors. 
Colors indicate the input temperature to the simulated 
SEDs: black is -- lOK and Hght gray is - 35 K. 



ture d ifference is less than the errors given in lRodon et al.l 
(|2010 t). The median absolute difference in /3- value is 0.2 and 
for all apertures ex cept for two the diff erences are less than 
the errors given in iRodon et al.l (|2010[ ) . 

So why does our interpr etation of the — Td rela- 
tionship diff e r from that of lAnderson et al.l ()2010D and 
iRodon et al. l (2010)? Both of their results were based on 
15 apertures from one field while our result here is based 
on 129 apertures from eight fields; we believe our results 
are more robust. Additionally, their (3 — Td fits were depen- 
dant on results toward the "interior" of the bubbles, which 
we have excluded from the present analysis. Our think- 
ing has evolved on the issue of the bubble interiors and, 
given the range of dust temperatures likely present there 
(see Section [5.4^ . it seems prudent to exclude the interi- 
ors from any analysis involving single-temperature fits. We 
note that the tem peratures for the "interior" apertures in 
IRodon et al.l (|201Q) have significantly higher uncertainties 
than those of other apertures in the field. We conclude that 
the differences between the present work and the previous 
analyses stem partly from differences in data processing, 
but also from which apertures were considered and which 
70 /Ltm data were used. This demonstrates the caution re- 
quired for future work on the 13 — Td relationship. 



5.3. The Three-Dimensional Nature of Bubbles 

There has recently been some debate on whether the types 
of Hii regions studied here are two-dimensional objects, 
"rings," or three-dimensiona l objects, "bubbles." T h e J = 
3 H> 2 i^CO observations in 'Bea umont fc Wilhaiiisl (|2010[ ) 
show little emission toward the bubble interiors. The lack of 
detected CO emission in t he bubble interiors ind icates that 
"bubbles" are flat objects. lAnderson et al.l(|201lh . however, 
suggest that bipolar bubbles would be much more common 
if all bubbles were two-dimensional, an effect that is not 
observed. Furthermore, these authors point out that the 
recombination line widths of bubbles that are not bi-polar 
should be larger since the ionized flows toward and away 



from the observer are largely unimpeeded; such larger line 
widths arc not observed. 

The Herschel data shown here also support the idea 
that bubble H ii regions are three-dimensional structures. 
Toward the center of all bubbles in our sample there is 
emission at all -ffersc/ieZ wavelengths. While at wavelengths 
< 100 /im, such emission comes mainly from a warmer 
dust component (see below), the longer wavelength emis- 
sion indicates the presense of cold dust. If such bubbles 
were in fact rings, we would not expect to f i nd co ld dust 
toward their central regions. Anderson et al. (l2010l) found 
support for the hypothesis of' Beaumont fc Williamgj ()2010t) 
- namely little emission toward center of RCW120. This 
observation, however, relied on an earlier processing of the 
Herschel data that removed emission around bright PDRs 
and thus from the bubble interiors. The current process- 
ing using Scanamorphos better reconstructs the flux around 
and within the bubble interiors. 

What percentage of the total emission comes from the 
bubble interiors? Using aperture photometry, we answer 
this question by comparing the flux from the bubble inte- 
riors against that of the "Entire" class defined earlier. The 
four angularly largest regions in our sample, Sh 104, W5-E, 
RCW79, and RCW120 permit the most accurate measure- 
ments, and we restrict the current analysis to these regions. 
We use the "Entire" apertures from earlier and define new 
apertures at each wavelength for the interior regions that 
trace the inner PDR boundary. Because the choice of back- 
ground aperture may affect our measurements, we define 
four background apertures for each wavelength of the four 
Hii regions. We perform the aperture photometry as before 
to determine the percentage of the total emission detected 
towards the bubble "interiors." 

Roughly 20% of the total emission from the bubbles 
in our sample comes from the interior. This indicates that 
toward the interior of the bubbles, we are detecting cold 
dust emission from the "front" and "back" sides of the 
bubble. The average percentage of emission from the in- 
terior for the four regions studied here at 100 /im, 160 /im, 
250 ^m, 350 ^m and 500 ^m is 26±0.7%, 20±6%, 17±7%, 
14 ± 8%, and 12 ± 7%, respectively. Much of the deviation 
in the above numbers comes from Sh 104, which has 10% 
more fiux from the interior at each wavelength compared 
to the averages. The exact interpretation of these numbers 
requires modeling that is beyond the scope of this paper. 
It is, nevertheless, suggestive that some bubble Hii regions 
are three-dimensional objects. 

A higher percentage of the emission at lower wave- 
lengths comes from the bubble interiors. Restricting our 
analysis to RCW79 and RCW120, for which we have 
Herschel Hi-Gal observations at 70 /.tm and MIPSGAL ob- 
servations at 24 /xm, we find that the contribution to the 
total flux from the interior is 60% at 24 ^m for RCW 120 
and 55% for RCW 79, and 33% at 70 /im for RCW 120 and 
69% for R CW 79. These 24 urn va lues are slightly higher 
than that of lDeharveng et alj (|2010( ). who found that about 
half the emission at 24 /j,m comes from the bubble interiors. 

5.4. The Warm Dust Component 

Nearly all regions in our sample have a "warm" dust 
component not accounted for in the current treatment 
(see Figure [7]) . This component has significantly less flux 
compared to the cold component. For a given flux, the 
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mass decreases with increasing temperature. The lower 
flux and higher temperature together imply that the warm 
component contains very little mass compared to the 
cold component. The warm compone nt emits strongly at 
24 ^m for bubb l e Hi i regions (see IWatson et al.l l2008t 
iDeharveng et al.l l201Clt) . but Herschel observations also 
show emission at longer wavelengths, especially 70 fim and 
100 /im, with the same spatial distribution seen at 24 /im 
(Figure [T51) . This indicates that some fraction of the emis- 
sion in the Herschel PACS bands is due to the warm 
component. This warm component has generally been at- 
tributed to s tochastically heated very small grains (VSGs; 
iDrainefc Lili20071 . 

Using the dust model DustEM (jComoiegne et al.ir201ltl , 
ICompiegne et al.l ()2010D show the SED Herschel Hi-Gal 
field at £ = 59°, a "diffuse" region and a "dense" region. 
In the "diffuse" field they found that the VSG component 
contributes on average 35% of the flux at 70 ^m, 13% of 
the flux at 100 /im, and 5% of the flux at 160 /im. In the 
"dense" field, these percentages decreased to 12%, 4%, and 
2%, respectively. Neither region contains active star forma- 
tion. Since nearly all of our measurements are toward dense 
structures, it seems likely that the latter numbers are more 
representative of our data. 

What effect does including data at 70 /im have on our 
derived temperatures ? Of the data used here to derive dust 
temperatures, the 70 /im data are most heavily affected by 
a possible warmer dust contribution, caused in part by the 
emission from VSGs. To answer this question, we perform 
two tests. In the first test, we recompute the temperature 
map for RCW 120 including the 70 /im data, and compare it 
to the map computed without the 70 /tm data. We find that 
including or excluding the 70 /im data point has minimal 
impact on the derived temperatures: the average change 
per temperature map pixel is ~ 4%, or 1.6 K on average. In 
our second test, we re-fit the fluxes extracted with aperture 
photometry for RCW 120. The temperatures change by ^ 
6% on average for RCW 120 apertures, compared to the 
temperatures derived when the 70 /im data are included 
when /3 is held to a value of 2.0. When /3 is allowed to vary, 
the derived temperatures change by ~ 4% and the derived /3 
values change by '^ 8%. In both cases, excluding the 70 /im 
data has the effect of decreasing the derived temperature 
by < 1 K on average, although the derived /3- values are on 
average neither higher nor lower. Excluding data at 70 /im, 
however, increases the fit uncertainties (cf. Figure [TT|) . We 
conclude that including data at 70 /im data has minimal 
impact on the derived temperatures and beta-values. 

We stress that the impact of excluding the 70 /im is less- 
ened because we also have data at 100 /tm. We re-fit the 
fiuxes extrated using aperture photometry for RCW 120 ex- 
cluding data at both 70 /im and 100 /tm. When /3 is held 
fixed to a value of 2.0, we find temperatures ^ 40% differ- 
ent than when the 70 /im and 100 /im data were included. 
For apertures with previously-derived temperatures greater 
than 30 K, the temperature difference is up to 100%. Many 
regions of the Hi-Gal survey will not have access to Herschel 
data at 100 /im. We caution that the derived temperatures 
for regions > 30K will not be well-constrained if 70 /im 
data are excluded and 100 /im data arc not available. 

What effect may the warmer dust component have on 
our derived temperatures? T^o answer this question, we re-fit 
the data for RCW 120 after decreasing the fiux to account 
for the warmer dust. To account for the largest possible con- 



tribution, we adjusted the fl uxes downward by the a verage 
values for the diffuse field in'Compiegn e et aTl (|2010f ) (35%, 
13%, and 5% for 70 /im, 100 /im, and 160 /tm) and re-fit 
the SED. We find that the derived temperatures change by 
~ 5% both when j3 is fixed and when it is allowed to vary, 
and that the derived /3-values change by ^ 5%. As men- 
tioned above, these probably over-estimate any contribu- 
tion from a warmer dust component for the dense structures 
here studied, and therefore the 5% temperature change 
should be regarded as an upper limit. We conclude that 
the warmer dust component, caused in part by VSGs, has 
minimal effect on the derived temperatures and j3-values. 



5.5. The Masses of Bubble HII Regions 

The collect and collapse process requires a massive shell 
of material for the formation of second-generation stars. 
With the calculated masses, we may test whether the de- 
rived mass is similar to that expected based on the size of 
the region. We show in Figure [TJ] the relationship between 
the total mass of the H ii region from the "Entire" aper- 
tures and the expected mass based on the physical size of 
the region, for a range of ambient densities from 100 to 
1, OOOcm^-^. We calculated the expected mass by assuming 
that all material overtaken by the ionization front during 
the expansion of the H ii region is still associated with the 
region, by neglecting the mass from the ionized plasma, and 
by assuming that the Hii regions are spherical. The ex- 
pected mass is therefore the dust density times the volume 
of the sphere times a gas-to-dust ratio of 100 (the bubble 
diameters are given in Table 1). The assumption that the 
bubbles are spherically symmetric introduces uncertainty 
to the calculations. It is nevertheless a good first-order as- 
sumption and is supported by our finding that the bubbles 
are likely three-dimensional (see Section [573)1 . 

Figure [U] shows that for most regions, the computed 
mass is similar to the expected mass for reasonable values 
of the ambient density. The largest exception is for W5- 
E, which has significantly less mass than expected. The 
expansion of W5-E may not have been restricted, as is im- 
plied by its open geometry. This region may have formed 
in a filament oriented east-west, and its expansion to the 
north and south may have progressed largely unimpeded 
(jDeharveng et al.ll2012t ). The opposite may be the case for 
RCW 71 and G332, which have more mass than expected 
and appear to have nearly complete PDRs. 



6. Conclusions 

We have analyzed the dust properties of eight Galactic H ii 
regions with a bubble morphology using Herschel data from 
the HOBYS and "Evolution of Interstellar Dust" key pro- 
grams. We used aperture photometry to determine the dust 
temperature and dust emissivity index f3 for regions of in- 
terest within the fields of these eight regions. We also com- 
puted maps of the dust temperature and the total column 
density. Most locations show contributions from two dust 
components. Here we only determine temperatures for the 
colder of these two components, which contains nearly all 
the dust mass. 

Our conclusions are: 

— All H II regions in the sample are well-characterized with 
mean dust temperatures of 26 K. 
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Fig. 13: The warm dust component in RCW79. The warm dust is prominent at 24 fim and decreases in intensity relative 
to the emission from the PDR at longer wavelengths. The size of all fields is 20' square. 
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Fig. 14: The expected swept up mass versus the calcu- 
lated mass for the eight bubble regions. The expected 
mass is computed for initial ambient densities from 100 to 
l,000cm~'^, assuming all material overtaken by the ioniza- 
tion front is now contained in the PDR of the region. The 
dashed line shows a one-to-one relationship. 



— The PDRs of the H ii regions have mean temperatures of 
26 K. The temperature is an important parameter for in- 
vestigating subsequent star formation since higher tem- 
peratures lead to higher values of the Jeans mass, all else 
being equal. At the wavelengths sampled by Herschel, 
bubble H ii regions are dominated by the emission from 
their PDRs - approximately 80% of the total flux is 
contained within the PDR. The remaining 20% of the 
flux is from the interior of the bubbles, which suggests 
that the bubbles are three-dimensional structures rather 
than rings. 

— The apertures can be divided into two statistically dis- 
tinct groups: a "warm" group averaging 26 K and a cold 
group averaging 17 K. The "warm" group contains aper- 
tures from the PDRs of H ii regions, as well as the to- 
tal emission from the H ii regions themselves while the 
"cold" group contains apertures from cold filaments. 

— The coldest regions of the fields are spatially coincident 
with infrared dark clouds (IRDCs); these have mean 
dust temperatures of 15 K. One third of the IRDCs here 
detected by Herschel have 160 /xm fluxes greater than 
their 100 /im fluxes. This criterion can be used to locate 
clouds with similar temperatures that, because of the 
lack of mid-infrared background, do not appear infrared 
dark. Filaments that are not detect in absorption at 
8.0 /im have mean dust temperatures of 19K. 



We estimate that not accounting for the warmer dust 
component and including the data at 70 /im affects the 
derived temperatures by at most 5% for the objects 
studied here. We stress that the inclusion of data at 
70 /im has less of an impact because we also have data 
at 100 /im - the exclusion of data at both 70 /im and 
100 /im results in large differences in temperature of 
^ 40%. This will be the situation for the majority of the 
Herschel Hi-Gal observations. In this case, the inclusion 
of data at 70 /xm is preferable to its exclusion. Color 
corrections employed here change the derived tempera- 
tures by ^ 3%. 

The total mass associated with these regions spans a 
wide range, from ~ 300 to ~ 10, 000 Mq. 
We find only weak support for a relationship between 
f3 and T^, in contrast to what has been found by pre- 
vious authors. The slight correlation between f3 and Td 
that we found is statistically similar to that caused by 
calibration and photometric uncertainties. This result 
differs from what was suggested by earlier analyses of 
the data shown here, reduced using a different process- 
ing method. That the P — Td relationship is so sensitive 
to the data processing method shows the difficulty in 
trying to understand it using Herschel data alone. 
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Table 3. Dust properties derived from aperture photometry 
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Shl04 


Entire 


74.765 


+0.624 


75.0 


24.4 + 0.8 


4200 


27.6 + 3.9 


1.7 + 0.3 


3300 






HII (UCHII) 


74.790 


+0.572 


2.7 


28.1 ±1.0 
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28.7 + 4.2 


2.0 + 0.3 


420 


I 




CS (CS) 


74.721 
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0.5 


24.1+0.7 


89 


25.0 + 3.2 


1.9 + 0.3 


82 






E. PDR 
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+0.565 


1.2 


24.9 + 0.8 


140 


24.4 + 3.0 


2.1 + 0.3 


150 


F 




N. PDR 1 (PDR 2) 


74.812 


+0.615 


0.9 


22.2 + 0.6 


120 


23.6 + 2.8 


1.8 + 0.3 


110 


F 




N. PDR 2 (PDR 3) 


74.818 


+0.657 


2.0 


24.0 + 0.7 


110 


24.3 + 3.0 


2.0 + 0.3 


110 


F 




W. PDR 1 (PDR 4) 


74.787 


+0.679 


2.6 


21.3 + 0.6 


300 


23.4 + 2.7 


1.7 + 0.3 


240 


F 




W. PDR 2 (PDR 6) 


74.728 


+0.671 


2.4 


23.8 + 0.7 


160 


24.7 + 3.1 


1.9 + 0.3 


150 


F 
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+0.624 


1.4 


25.0 + 0.8 


160 
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F 
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1.2 
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49 
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50 


F 
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+0.705 
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31 
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N. Filament 2 
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18.5 + 1.7 
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N. Filament 3 
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+0.618 


1.3 


17.9 + 0.4 


110 
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1.7 + 0.3 


89 






N. Filament 4* (PDR 1) 


74.836 


+0.583 
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19.9 + 0.5 
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1.7 + 0.3 


140 
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I 
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23.8 + 0.7 
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I 
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2.4 


21.0 + 0.6 


43 


22.6 + 2.5 


1.8 + 0.3 


37 
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+1.560 


4.4 


23.5 + 0.7 


440 


25.6 + 3.3 


1.8 + 0.3 


370 






N. PS 1 (BRC13) 


138.142 


+1.684 


1.6 


22.4 + 0.6 


60 


24.2 + 2.9 


1.8 + 0.3 


51 






W. PS 1 


137.793 


+1.461 


1.7 


23.4 + 0.7 


49 


25.0 + 3.1 


1.8 + 0.3 


43 






W. PS 2 (BRC12) 


137.505 


+1.397 


1.3 


21.9 + 0.6 


18 


22.6 + 2.6 


1.9 + 0.3 


17 
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137.537 
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1.5 


21.2 + 0.6 


40 
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1.7 + 0.3 


32 
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25.8 + 3.4 


1.6 + 0.3 


580 


F 




E. PDR 2 


138.319 


+1.671 


3.7 


21.9 + 0.6 


54 


23.7 + 2.8 


1.8 + 0.3 


46 


F 




N. PDR 1 (BRC13) 


138.187 


+1.727 


24.7 


21.6 + 0.6 


400 


23.0 + 2.6 


1.8 + 0.3 


350 


F 




N. PDR 2 


138.031 


+1.762 


3.9 


22.5 + 0.6 


27 


21.0 + 2.2 


2.2 ± 0.3 


32 


F 




N. PDR 3 


137.921 


+1.664 


19.2 


21.2 + 0.6 


170 


24.3 ± 2.9 


1.6 + 0.3 


130 


F 




N. PDR 4 


137.897 


+1.920 


43.2 


21.7 + 0.6 


160 


20.1 + 2.0 


2.3 + 0.3 


200 


F 




W. PDR 1 


137.930 


+1.119 


196.1 


23.1+0.7 


1100 


24.3 + 3.0 


1.9 + 0.3 


1000 


F 




W. PDR 2 (BRC12) 


137.553 


+1.344 


6.9 


23.8 + 0.7 


67 


24.3 + 3.0 


1.9 + 0.3 


64 


F 




N. Filament 1 


138.256 


+1.850 


1.1 


18.4 + 0.4 


7.6 


19.4 + 2.0 


1.8 + 0.4 


6.6 






N. Filament 2 


137.738 


+1.919 


5.4 


15.8 + 0.3 


54 


17.0 + 1.4 


1.7 + 0.4 


45 






E. Filament 1 


138.596 


+1.707 


12.4 


17.5 + 0.4 


67 


19.1 + 1.8 


1.7 + 0.3 


54 






E. Filament 2 


138.594 


+1.584 


4.7 


18.3 + 0.4 


79 


18.7 + 1.7 


1.9 + 0.3 


75 




Sh241 


Entire 


180.886 


+4.114 


16.5 


22.9 + 0.7 


1800 


27.2 + 3.8 


1.5 + 0.3 


1300 






HII 


180.828 


+3.993 


2.6 


24.8 + 0.8 


350 


27.6 + 3.9 


1.7 + 0.3 


290 


I 




N. PS* 


180.820 


+4.189 


0.4 


19.2 + 0.5 


32 


21.0 + 2.2 


1.7 + 0.3 


26 






W. PS* 


180.859 


+4.032 


0.5 


16.6 + 0.4 


93 


20.0 + 1.9 


1.3 + 0.3 


58 






S. PS* 


180.918 


+4.017 


0.9 


18.4 + 0.4 


110 


19.5 + 1.9 


1.8 + 0.3 


98 






N. PDR 1 


180.856 


+4.122 


0.8 


23.8 + 0.7 


150 


25.8 + 3.4 


1.8 + 0.3 


130 


F 




N. PDR 2 


180.858 


+4.106 


0.7 


24.1+0.8 


150 


28.1+4.0 


1.6 + 0.3 


110 


F 




W. PDR 


180.873 


+4.092 


1.3 


22.5 + 0.7 


560 


25.8 + 3.3 


1.6 + 0.3 


420 


F 




W. Filament 1* 


180.855 


+4.069 


0.4 


19.2 + 0.5 


39 


23.9 + 3.3 


1.3 + 0.4 


24 






W. Filament 2* 


180.840 


+4.044 


5.4 


18.5 + 0.4 


540 


20.0 + 2.0 


1.7 + 0.3 


450 






W. Filament 3* 


180.866 


+4.053 


0.3 


15.4 + 0.3 


110 


17.6 + 1.6 


1.5 + 0.4 


79 






W. Filament 4* 


180.874 


+4.047 


0.7 


16.1+0.4 


140 


19.7 + 1.9 


1.2 + 0.3 


82 






S. Filament* 


180.915 


+4.040 


0.9 


17.3 + 0.4 


100 


18.4 + 1.9 


1.8 + 0.4 


87 




RCW71 


Entire 


302.799 


+1.287 


21.6 


30.3 + 1.2 


240 


30.8 + 5.0 


2.0 + 0.3 


230 






S. PS 1* 


302.842 


+ 1.225 


0.8 


15.3 + 0.3 


37 


15.8 + 1.4 


1.9 + 0.4 


34 






S. PS 2* 


302.809 


+1.233 


0.9 


18.2 + 0.4 


20 


19.2 + 2.0 


1.8 + 0.4 


18 






N. PDR 


302.802 


+1.309 


1.6 


34.5 + 1.6 


13 


30.2 + 5.0 


2.3 + 0.4 


16 


F 




W. PDR 


302.780 


+ 1.284 


0.7 


30.7 + 1.2 


15 


30.2 + 4.9 


2.0 + 0.3 


15 


F 




S. PDR 1 


302.784 


+1.259 


1.1 


28.3 + 1.0 


20 


27.5 + 4.0 


2.1 + 0.3 


21 


F 




S. PDR 2 


302.807 


+1.261 


0.8 


26.4 + 0.9 


24 


28.2 + 4.1 


1.8 + 0.3 


21 


F 




E. PDR 1 


302.822 


+1.291 


0.6 


34.9 + 1.7 


8.7 


33.4 + 6.2 


2.1 + 0.4 


9.3 


F 




E. PDR 2 


302.810 


+1.277 


0.6 


33.2 + 1.5 


19 


30.4 + 4.9 


2.2 + 0.3 


22 


F 




E. PDR 3 


302.823 


+1.273 


0.5 


32.6 + 1.4 


6.9 


31.6 + 5.6 


2.1 + 0.4 


7.3 


F 




W. Filament 1* 


302.718 


+1.319 


2.3 


16.1+0.4 


40 


16.9 + 1.6 


1.8 + 0.4 


35 






W. Filament 2* 


302.732 


+1.307 


0.5 


18.3 + 0.5 


5.8 


18.2 + 2.5 


2.0 + 0.6 


6.0 
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W. Filament 3* 


302.720 


+1.271 


0.6 


14.3 + 0.3 


12 


14.3 + 1.7 


2.0 + 0.6 


12 




W. Filament 4* 


302.734 


+1.276 


0.5 


17.3 + 0.4 


8.9 


16.9 + 2.0 


2.1 + 0.5 


9.3 




N. Filament* 


302.859 


+1.388 


0.8 


12.8 + 0.6 


9.8 


9.6 + 1.9 


3.9 + 1.6 


29 




RCW79 Entire 


308.709 


+0.628 


182.4 


26.9 + 0.9 


11000 


28.8 + 3.8 


1.8 + 0.3 


9700 




HII 


308.750 


+0.548 


6.8 


27.5 + 0.9 


1700 


25.0 + 2.0 


2.2 + 0.2 


2100 


I 


W. PS 1 (Cond. 4) 


308.645 


+0.648 


1.4 


25.1+0.8 


720 


26.0 + 2.4 


1.9 + 0.2 


660 




W. PS 2 (Cond. 6) 


308.673 


+0.637 


1.1 


29.7 + 1.1 


73 


27.7 + 3.1 


2.2 + 0.2 


84 




N. PDR 1 (Cond. 5) 


308.729 


+0.725 


2.6 


23.4 + 0.6 


440 


21.8 + 1.7 


2.2 + 0.2 


520 


F 


N. PDR 2 


308.793 


+0.648 


14.0 


29.2 + 1.0 


660 


21.1 + 1.6 


2.8 + 0.2 


1300 


F 


W. PDR 1 


308.630 


+0.558 


3.6 


34.3 + 1.5 


170 


27.8 + 2.9 


2.4 + 0.2 


250 


F 


W. PDR 2 


308.613 


+0.607 


4.2 


36.1 + 1.6 


250 


22.6 + 1.8 


2.9 + 0.2 


610 


F 


W. PDR 3 (Cond. 7) 


308.673 


+0.620 


0.9 


27.1 + 1.0 


69 


28.2 + 3.2 


1.9 + 0.2 


64 


F 


S. PDR (Cond. 3) 


308.690 


+0.533 


6.6 


25.2 + 0.8 


1400 


25.8 + 2.3 


1.9 + 0.2 


1300 


F 


2nd PDR E. 


308.708 


+0.470 


23.5 


24.3 + 0.7 


2000 


24.4 + 3.0 


2.0 + 0.3 


2000 


F 


2nd PDR W. 


308.522 


+0.585 


14.7 


23.1+0.7 


610 


20.5 + 2.1 


2.4 + 0.4 


800 


F 


N. Filament (Cond. 8) 


308.725 


+0.774 


6.9 


21.7 + 0.5 


950 


18.3 + 1.1 


2.5 + 0.2 


1400 




E. Filament* 


308.892 


+0.698 


7.9 


24.0 + 0.7 


130 


22.1+2.4 


2.2 + 0.3 


160 




RCW82 Entire 


310.980 


+0.408 


62.2 


24.6 + 0.7 


3400 


21.3 + 2.2 


2.4 + 0.3 


4700 




E. PS 


311.034 


+0.382 


1.8 


21.1+0.5 


430 


20.4 + 1.4 


2.1 + 0.2 


470 




N. PDR 1 


311.003 


+0.462 


2.2 


25.3 + 0.8 


200 


22.7 + 1.9 


2.3 + 0.2 


250 


F 


W. PDR 1 


310.948 


+0.424 


2.8 


29.7 + 1.0 


130 


20.4 + 1.9 


3.0 + 0.3 


280 


F 


W. PDR 2 


310.928 


+0.445 


2.1 


25.8 + 0.8 


190 


20.5 + 1.7 


2.6 + 0.3 


310 


F 


S. PDR 


310.987 


+0.359 


4.8 


25.5 + 0.7 


340 


19.8 + 1.4 


2.7 + 0.2 


590 


F 


E. PDR 1 


311.026 


+0.413 


1.5 


25.3 + 0.7 


160 


19.9 + 1.4 


2.6 + 0.2 


270 


F 


E. PDR 2 


311.013 


+0.430 


1.1 


23.0 + 0.6 


110 


19.3 + 1.4 


2.5 + 0.2 


160 


F 


N. Filament* 


311.060 


+0.516 


2.6 


21.9 + 0.6 


79 


19.0 + 1.8 


2.5 + 0.4 


110 


F 


W. Filament* 


310.850 


+0.486 


4.0 


18.4 + 0.4 


260 


17.5 + 1.4 


2.2 + 0.3 


300 




S. Filament 1* 


310.927 


+0.277 


2.7 


18.0 + 0.6 


57 


17.6 + 2.3 


2.1 + 0.5 


60 




S. Filament 2* 


310.896 


+0.379 


2.5 


19.3 + 0.5 


97 


16.5 + 1.6 


2.6 + 0.4 


140 


F 


E. Filament 1* 


311.086 


+0.429 


5.5 


19.7 + 0.5 


300 


19.4 + 1.5 


2.0 + 0.3 


310 




E. Filament 2* 


311.052 


+0.394 


1.1 


13.0 + 0.3 


290 


12.1 +0.9 


2.4 + 0.4 


380 


F 


G332.5-0.1 Entire 


332.539 


-0.141 


13.6 


25.9 + 0.8 


4500 


26.1+3.3 


2.0 + 0.3 


4500 




HII 1 


332.544 


-0.126 


0.9 


29.2 + 1.0 


590 


29.1+3.1 


2.0 + 0.2 


590 


1 


HII 2 


332.445 


-0.178 


4.4 


24.0 + 0.7 


330 


19.3 + 1.4 


2.6 + 0.2 


540 


1 


W. PS 


332.603 


-0.167 


1.6 


20.0 + 0.5 


1200 


20.4 + 1.4 


1.9 + 0.2 


1100 




E. PS 


332.467 


-0.130 


1.8 


20.7 + 0.5 


850 


18.5 + 1.2 


2.4 + 0.2 


1100 




E. PDR 


332.560 


-0.149 


1.1 


22.6 + 0.6 


650 


20.9 + 1.5 


2.2 + 0.2 


770 


F 


S. PDR 


332.526 


-0.155 


1.4 


25.5 + 0.8 


650 


23.7 + 2.0 


2.2 + 0.2 


750 


F 


E. Filament* 


332.648 


-0.188 


3.2 


15.9 + 0.3 


700 


15.1 + 1.1 


2.2 + 0.3 


810 


F 


S. Filament* 


332.404 


-0.150 


5.0 


19.6 + 0.5 


820 


21.7 + 1.9 


1.7 + 0.3 


660 


F 


W. Filament 


332.501 


-0.119 


1.5 


22.5 + 0.6 


710 


19.9 + 1.3 


2.3 + 0.2 


930 


F 


RCW 120 Entire 


348.251 


+0.481 


100.3 


23.9 + 0.6 


2100 


26.0 + 2.4 


1.7 + 0.3 


1800 




Cond. 1 PS 


348.184 


+0.482 


0.7 


20.9 + 0.5 


300 


22.6 + 1.6 


1.8 + 0.2 


260 




Cond. 1 (Cond. 1) 


348.180 


+0.480 


4.7 


21.0 + 0.5 


760 


23.5 + 1.6 


1.7 + 0.2 


600 


F 


S. PDR 1 (Cond. 2) 


348.226 


+0.413 


2.5 


24.1+0.6 


160 


25.5 + 2.0 


1.8 + 0.2 


140 


F 


E. PDR 1 (Cond. 3) 


348.272 


+0.403 


2.6 


23.4 + 0.6 


88 


23.4 + 1.9 


2.0 + 0.2 


88 


F 


E. PDR 2 (Cond. 4) 


348.302 


+0.429 


3.0 


22.6 + 0.6 


110 


22.3 + 1.7 


2.0 + 0.2 


110 


F 


W. PDR 1 (Cond. 7) 


348.199 


+0.531 


0.7 


24.8 + 0.7 


20 


23.4 + 1.9 


2.1 + 0.2 


23 


F 


S. PDR 1 (Cond. 10) 


348.203 


+0.426 


1.0 


26.2 + 0.8 


38 


25.8 + 2.2 


2.0 + 0.2 


39 


F 


W. PDR 2 (PDR West) 


348.266 


+0.556 


0.8 


31.6 + 1.2 


6.0 


20.7 + 1.6 


3.0 + 0.2 


14 


F 


E. PDR 3 (PDR East) 


348.329 


+0.464 


1.2 


25.3 + 0.7 


31 


20.8 + 1.4 


2.5 + 0.2 


47 


F 


W. IRDC 1* 


348.179 


+0.600 


1.5 


12.7 + 0.3 


30 


12.0 + 0.9 


2.3 + 0.4 


37 


F 


W. IRDC 2* 


348.165 


+0.662 


2.3 


14.0 + 0.3 


30 


13.1+0.8 


2.3 + 0.3 


37 


F 


W. IRDC 3* 


348.196 


+0.768 


1.8 


13.5 + 0.2 


64 


14.1 + 1.1 


1.8 + 0.4 


55 


F 


W. IRDC 4* 


348.234 


+0.591 


1.8 


17.6 + 0.4 


28 


19.6 + 2.0 


1.6 + 0.4 


22 


F 


N. IRDC 1* (Cond. 6) 


348.292 


+0.644 


2.0 


15.0 + 0.3 


160 


15.4 + 0.8 


1.9 + 0.2 


150 


F 


N. IRDC 2* (Cond. 9) 


348.345 


+0.652 


1.5 


13.4 + 0.2 


110 


12.5 + 0.6 


2.3 + 0.2 


140 


F 


N. IRDC 3* 


348.297 


+0.682 


1.4 


13.1+0.6 


26 


11.9 + 1.5 


2.3 + 0.4 


36 


F 


N. IRDC 4* 


348.305 


+0.703 


0.8 


12.2 + 0.2 


38 


9.9 + 0.4 


3.2 + 0.3 


88 


F 


N. IRDC 5* 


348.325 


+0.716 


1.4 


11.5 + 0.4 


73 


10.1 + 1.1 


2.5 + 0.4 


120 


F 


N. IRDC 6* 


348.476 


+0.583 


1.3 


11.8 + 0.5 


61 


11.6 + 1.4 


2.1 + 0.4 


65 


F 


E. IRDC* (Cond. 5) 


348.375 


+0.487 


10.6 


15.9 + 0.3 


520 


16.3 + 0.9 


1.9 + 0.2 


490 


F 


S. Filament 1 (S. Filament 1)* 


348.121 


+0.473 


5.4 


16.6 + 0.3 


160 


16.6 + 1.0 


2.0 + 0.2 


160 




S. Filament 2 (S. Filament 2)* 


348.134 


+0.512 


4.4 


15.5 + 0.3 


140 


15.2 + 0.9 


2.1 + 0.3 


140 




S. Filament 3 (Cond. 8)* 


348.149 


+0.421 


3.3 


16.6 + 0.3 


160 


16.1+0.9 


2.1 + 0.2 


180 
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'^ The p arenthetical names are from previous wo rks: |Rod6n et alj (120101 ) fo r Sh 104. [ Sharplessl (Il959l ) and lSugitani et aU (llQQlh 
for W5-E. IZavagno et all ((200i) for RCW79, and Zavasno ct al. (2007) and lAnderson et all (|2010l ') for ROW 120. The names of 
apertures where we used a local background value arc marked with an asterisk. 

'' The classifications are: "E" for entire, "Hll" for other Hll regions in the field, "PS" for point source, "PDR" for apertures along 
the photodissociation regions, "F" for filaments, and "F (ID)" for IRDC filaments seen in absorption at 8.0 /im. 
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Appendix A: Comments on Individual Regions 

Below we comment briefly on the eight bubble regions in 
our sample. The masses and temperatures quoted are those 
of the /3 = 2.0 trial. 



A.l. Shl04 

We find a total associated mass for Sh 104 of ^ 4, 000 Mq, 
which is in the middle of the range for the H ii regions inves- 
ti gated here. This calculat ed mass is two-thirds that found 



Deharveng et al.l ([2003) using molecular gas tracers. As 



mentioned previously, the interior of Sh 104 has more emis- 
sion at all wavelengths than the other regions in the sam- 
ple. Roughly 30% of the total FIR emission from the region 
comes from the interior region. We hypothesize that this is 
a result of Sh 104 being nearly complete - some of the ra- 
diation that would otherwise escape is trapped within the 
bubble. 

In Figure[TJ there is a clear difference in temperature be- 
tween Sh 104 itself and the local ambient medium. We find 
that its PDR is well-characterized by temperatures of 25 K. 
Local filaments are for the most part ^ 20 K, although "N. 
Filament 1" and "N. Filament 5" appear to be significantly 
colder at 15 K and 17 K respectively. The temperature map 
of Sh 104 (FigureHl) is relatively smooth but there is a trend 
for hotter temperatures toward the southwest, and colder 
temperatures toward the north-west. The warmest region 
in the field is the separate Hii region at {£, b) = (73.790, 
-1-0.572). It has a mean dust temperature of 28 K and a to- 
tal associated mass of ~ 400 Mq. The other source found 
along the PDR, "CS" has a mass of ~ 100 Mq. 

The column density map in Figure [S] shows the highest 
values, up to lO^^cm"^, along the PDR. Local filaments 
have mean column densities of > lO^^cm"^. There are a 
number of prominent filaments that begin on the border of 
Sh 104 and lead radially away; these are prominent in the 
column density map. 



A.2. W5-E 

The most surprising aspect of W5-E found here is that 
its dust temperatures vary little across the field. Figure [1] 
shows a narrow range of Herschel colors, which leads 
to relatively small variations in dust temperature. The 
dust temperature of the entire region is 24 K. The cold- 
est regions of the field are to the east ("E. Filament 1" 
and "E. Filament 2") and north ("N. Filament 1" and 
"N. Filament 2") - they are between 16 and 18 K. There 
are other cold locations along the PDR, most notably be- 
hind the ionization front in "N. PDR 1" (BRC13) and 
"E. PDR 1 (BRC14)." The total associated mass of W5-E, 
~ 8, 000 Mq, places it as the second most massive in the 
sample; it is the largest region in the sample in terms of 
physical diameter. 

There are three other Hii regions in the field of W5- 
E. With the exception of the separate Hii region on the 
border of in RCW79, these are the most massive secondary 
Hii regions in our sample. The bipolar Hii region Sh201 
at {£, b) = (138.481, -hi. 637) is the hottest location in the 
field at 27 K and stands out clearly in Figure 2] Its total 
associated mass is ~ 800 Mq. The other Hii regions on 
the eastern border are ■^ 25K; they have total associated 
masses of ^ 600 Mq. 



W5-E is unique in our sample in that it displays a 
number of features obviously affected by the impinging 
stellar radiation. These features are discussed in detail in 
iDeharveng et al.l (|2012D . The bright-rimmed clouds BRC12 
and BRC13, which harbor embedded sources, have mean 
temperatures of 24 K and 22 K, respectively, while the dust 
associated with the embedded sources has a temperature 
of ^ 22 K. The embedded sources are conspicious in the 
temperature map (Figure S]) as higher temperature regions 
within cold temperature clumps. The pillars detected with 
Spitzer seen to the southwest are also cool and have an aver- 
age temperature of ^ 20 K. Many of these pillars have em- 
bedded stars forming at their "tips" detected by Herschel. 
This is consistent with a scenario where the majority of the 
pillar is shielded from impinging ratiation and is thus able 
to maintain the cold temperatures necessary for subsequent 
star formation. 



A.3. Sh241 

The temperature map of Sh 241 in Figure S] shows relatively 
little variation across the field. The large high-temperature 
region to the north is probably not real as it is a region of 
very low intensity emission. The bubble H ii region itself is 
23 K. The prominant separate compact H ii region is slightly 
warmer, 25 K while filaments average ~ 17 K. The total 
mass for Sh241 is ^ 2, 000 Mq while that of the compact 
Hii region is ~ 300 Mq. Sh241 is ringed by a region of cold 
emission to the north of temperature ~ 18 K. This region 
has a higher than average column density, > lO^^'^cm"^. 
The colder emission appears to be part of a large filament 
running East- West that also contains another Hii region 
and what appears to be cold protostars. 



A.4. RCW71 

As a whole, RCW71 is the warmest region in our sample at 
30 K. RCW 71 is also the least massive region in the sample, 
with a total associated mass of just ^ 200 Mq . This low 
mass estimate may be due to an inaccurate distance. As 
mentioned in SQl the kinematic and spectroscopic distances 
for this source do not agree. We note, however, that the 
PDR of RCW 71 has only marginally higher column density 
values than the background (Figure ^ , and therefore the 
low associated mass estimate may be real. 

The temperature map of RCW 71 in Figure S] shows 
high temperatures in the PDR up to '^ 35 K (in apertures 
"E. PDR 1" and "E. PDR 2"). Material has accumulated 
into cold massive filaments most prominently to the east, 
but also to the west. The coldest of these, "N. Filament" 
has a temperature of 13 K - it is one of the coldest regions in 
the sample. The field of RCW 71 also has one of the coldest 
point sources, "S. PS 1." It has a temperature of 15 K. 



A.5. RCW 79 

RCW 79 has a total associated mass of ~ 10, 000 Mq, 
which makes it the most massive in our sample. This mass 
does not include the secondary PDRs to the south, which 
are themselves ~ 3, 000 Mq combin ed. Our mass es t imate s 
are considerably higher than that of lZavagno et al.l (|2006[ ). 
who found 2, 000 Mq for the entire region using 1.2 mm ob- 
servations. Their observations, however, may not have been 
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sensitive to the more diffuse emission detected by Herschel. 
The column density map in Figure[6]has values in the PDR 
of up to ^ 10^^-^ cm~^ and lo oks very similar to the 1.2 mm 
maps in lZavagno et al.l (J2006f) . Along the south-eastern bor- 
der there is another Hii region that has a total associated 
mass of ^ 2, 000 Af©. This is the most massive secondary 
H II region in our sample - it is significantly more mas- 
sive than RCW 71 and approximately the same mass as 
RCW120. 

RCW 79 is interesting in that there is significant patches 
of cold material ringing the PDR to the south and to the 
west. These colder areas are apparent in Figure SI The 
coldest objects in the field, the "N. Filament" and the "E. 
Filament" are 22 K and 24 K, respectively. There is ongo- 
ing star formation in the field. Within "E. Filament" there 
are two compact objects and there are compact sources de- 
tected in the "S. PDR" aperture. 

A.6. RCW 82 

Like RCW 71, the PDR of RCW 82 is filamentary, but its 
temperature of 25 K is cooler than that of RCW 71. The 
total associated mass of RCW 82 is - 3, 00 M©. This is 
signifi cantly less than what was found by iPomares et al.l 
(|2009l) from CO observations. The temperature in the field 
of RCW 82 seen in Figure E] shows little variation. To the 
west, RCW 82 has what appears to be a second PDR 
("W. PDR 2"). The IRDCs seen in the field range from 
~ 20K for "N. Filament" to < 15K for "E. Filemant 2". 
As for Sh 104, there are a number of cool filaments lead- 
ing radially away from RCW 82; these have column density 
values of - lO^^'-^ cm~2. The eastern PDR of RCW 82 has 
a number of cold locations. These are evident in the col- 
ors seen in Figure [2J and also in the temperature map in 
Figure HI 

A.7. G332.5-0.1 

G332.5— 0.1 is angularly small and located along a promi- 
nent IRDC. The total associated mass of G332.5— 0.1 is 
~ 4, 000 Mq which is on the higher end of the range of H ii 
regions masses studied here. The IRDC itself has a temper- 
ature of up to ^ 20 K, which is warmer than other IRDCs 
here studied. Indeed, as seen in Figure HI the IRDC is not 
well-separated in temperature from the surrounding region. 
The column density distribution (FigurelB]) shows the IRDC 
prominently - it has column density values of ~ 10^^ cm~^. 
Within the IRDC there are two prominent condensa- 
tions: "W. PS" and "E. PS". These condensations have 
masses of ~ 1, 000 Mq. The PDR has massive condensa- 
tions to the north, south, and east. The northern conden- 
sation contains an H ii region, with a total associated mass 
of of ~ 600 Mq. This location is the hottest in the field, 
29 K. 



A.8. RCW120 

The wide range of colors seen Figure [1] indicates that 
RCW 120 also has a large range of dust temperatures, com- 
pared to the other H ii regions in the sample. RCW 120 has 
the largest concentration of IRDCs in its surroundings of 
any object in our sample. These IRDCs are the coldest 
objects found in our analysis, averaging 13 K. The largest 



of the IRDCs, "E. IRDC," has a mass of - 500 Mq and 
has numerous protostars detected within it. This IRDC 
has a temperature derived from aperture photometry of 
16 K. This temperature is almost certainly affected by the 
numerous embedded protostars which have not been re- 
moved when calculating temperatures. The temperature 
map (Figure 2]) shows temperatures of < 15 K for quiescent 
parts of the cloud where there are no detected protostars. 
There are two large protostars detected in the "E. IRDC" 
aperture that are visible in the temperature map - they 
have temperatures of ^ 18 K. 

A point source is detected within C ondensation 1 at 
70 ^m and 100 /im (jZavaeno et al.ll2010i ) for which we find 
a total associated mass of ~ 300 Mq. The entire condensa- 
tion has a mass of ~ 800 Mq. It is easily the most mas- 
sive condensation in the field. In fact, it's mass is one 
third that of the total associated mass for the RCW 120 
region, which is only ^ 2, 000 M© . Using the same 870 /xn i 
APEX-LABOCA data used here, iDeharveng et al.l (|2009t) 
calculated a mass for the source within Condensation 1 of 
140-250 Mq and a mass for the entire condensation of 460- 
800 Mq ; the range of values comes from estimates of the 
dust temperature from 30 K to 20 K. 

The column density map of RCW 120 (Figure [5]) high- 
lights the numerous filamentary structures. These are 
prominent in the source, but are also detected leading away 
from the "E. IRDC" region. The PDR and the local IRDCs 
have roughly the same value of A^h, ~ 10^^ cm~^. 

lAnderson et al.l (|2010f ) found that for RCW 120, the ion- 
ization front appears "patchy" and that there are numer- 
ous locations where hot dust is detected beyond the ion- 
ization front. They hypothesized that these locations rep- 
resent holes in the PDR where radiation may escape and 
heat the ambient medium. Such radiation may cause an 
increase in pressure which would aid in collapsing existing 
pre-stellar clumps. We confirm their results here; the same 
three warmer locations can again be identified beyond the 
ionization front of RCW 120. It is interesting, however, that 
RCW 120 remains the clearest example of this phenomenon. 
This may be due to its proximity to the Sun. 
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